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LAYERS 
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Institute of Physics, Polish Academy of Science, Warsaw 


( Received September 26, 1955) 


The microstructure of PbS layers has been investigated by means of electron diffrac- 
tion and electron microscopy. The PbS layers were obtained by evaporating in vacuo. 
The specimens to be examined were prepared by the method of formvar pseudo-replicas 
or directly by evaporation on thin films of formvar, collodion of aluminium. In general, 
the PbS layers consist of crystals 200—300 А in size; an influence of the substratum 
on the size or shape of these crystals has not been observed. The diffraction patterns of 
layers deposited directly by evaporation correspond to a regular face-centred lattice 
of the NaCl type. The diffraction patterns of the pseudo-replicas show characteristic 
departures from this structure consisting of a decrease in the diameter of the 200 ring 
and an increase in the diameter of the 220 ring as compared with, the diameters of the 
remaining rings. 


We have investigated the structure of photoconductive PbS layers with the aid 
of an electron microscope and electron diffraction. These layers are formed in a hard 
glass envelope (Fig. 1) having the shape of a flat cylinder 35 mm in diameter and 15—20 
mm high with two tungsten electrodes Ë sealed in glass. On one of the walls of the enve- 
lope are painted graphite electrodes which extend to the electrodes E. The envelope is 
connected with the vacuum plant by means of tube H. Photocells from which specimens _ 
for the microscope have been prepared were further provided with cylinders about 15 ` 
mm in diameter and 50 mm long, serving for introducing the specimen holders 
(Fig. 2). The formation process of a layer consists in two evaporations: first of lead 
tiosulphate, which is poured into the envelope through the tube H before the latter . 


15 connected to the vacuum plant, while the layer is condensed on wall B; the second 


by heating wall B with a burner, this layer is re-evaporated on wall А, which is intensely 
cooled with an air stream. In order to intensify the cooling of this wall, a piece of glass 
tube C is soldered ‚to the envelope. | 

In the first series of experiments, the results of which were given in a previous 
work (Feltynowski et al. 1954), a polished glass plate 1 mm thick has been placed on 
the wall А, hiding half of this wall. When the photocell was made and its photoconduc- 
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tive properties controlled, the envelope was cut through and a 0.75% dioxan solution 
of formvar was poured over the plate covered with the PbS layer. The dioxan soon 
evaporated and the formvar film thus formed (about 100 À thick) was detached from the 
glass with the help of a paper band coated with a layer of glue. The PbS layer was 
detached together with the formvar film, forming a so-called pseudo-replica. After 
being immersed in warm water, the form- 
var film with the PbS layer rose to the 
| | surface and was taken out on a specimen 
holder provided with a grid. The speci- 
men thus prepared was examined by 
means of diffraction and with the aid 
of the microscope. 


Specimen 
holders 


Fig. 1 Fig. 2 


PbS has а face-centred cubic structure, yet the electron diffraction pattern ob- 
tained by us (Phot. 1) shows deviations from this structure. On the basis of the data in 
Table 1, where the diameters of the rings, the intensity of the latter and the plane in- 
dices ascribed to them on the assumption of face-centred cubic structure are given, 
we find that these deviations consist in: i 

(1) the shifting of the ring 200 towards the middle by 5.8%, 

(2) the shifting of the ring 220 outwards by 4%. 

Moreover, there appear several additional rings with considerably fainter intensities, 
from which rings No. 3,7 and 15 could be ascribed to PbS, as their position corresponds 
better with the typical PbS structure. The remaining faint rings derive from admixtures. 
This experiment was repeated three times with three different photocells. The PbS 
pattern from the pseudo-replicas always showed the same characteristic deviations. 

The specimens thus prepared were also examined under the microscope (Phot. 
2, Запа 4). It was found that in the layer there appear crystals which are not equally 
distributed, but tend to form clusters. Observations showed that the shapes and the 
sizes of these crystals were different for various photo-cells; for instance, in one (Phot. 2), 
the crystals were longitudinal, 900—3000 А long and 250—1000 À wide, in another 
(Phot. 3 and 4), however, they were considerably smaller (ca 400 À) and tended to 
form chains (Phot. 3). | 

In the second series of experiments the layer was re-evaporated directly on speci- 
men holders placed instead of the glass plate, and introduced there only after the 
first evaporation. The films on the specimen holders were of different kinds: 
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CABBIE 1 
No D Intensity Plane indices ve D 
mm hkl سس تا‎ 
1 8.50 very strong ЛЕТ 4.91 
2 9.23 very strong 200 4.62 
3 9.83 faint (2 0 0) (4.92) 
4 10.67 faint 
5 011 faint 
6 12.25 faint 
Т 13.39 faint (2 2 0) (4.74) 
8 14.40 very strong 220 5.10 
9 14.89 very faint 
10 16.37 strong gi 4.94 
11 16.38 very faint 
12 17.20 medium 28282 4.97 
13 17.94 very faint 
14 18.63 faint 400 4.66 
15 19.55 very faint (4 0 0) (4.89) 
16 19.97 very faint 
17 21.38 strong 891 . 4,90 
18 21.76 medium < 420 4.87 
19 22.06 very faint 
20 22.72 faint 
21 23.34 very faint 
22 23.98 medium 422 4.90 
23 24.48 faint 
24 25.15 rather marked 3 5 4.84 
i 1 | 
25 25.40 faint 
26 25.94 very faint 
27 26.58 faint 
28 21.52 medium 440 4.86 
29 28.36 very faint 
30 29.08 faint 531 4.92 
31 29.60 faint l 0 | 4.93 
442 
32 32.90 medium 622 4.96 
33 33.70 medium 444 4.86 


collodion, formvar and aluminium. The essential difficulty consisted in putting the 
specimen holders into the envelope after the first evaporation, without heating the 
condensed layer. For this purpose the envelope was provided with a large tube D and 
suitable glass rods fór maintaining the specimen holders inside the envelope. 

The diffraction patterns obtained were: 

(1) quite in accordance with the PbS structure of the NaCl type (Phot. 5, Table 2); 

(2) in accordance with the PbS structure but having additional rings (Phot. 6, 
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TABLE 2 


D 


Intensity Dress 


1 strong 

2 very strong 

3 very strong 5.45 
4 18.0 strong 5.43 
5 18.8 rather marked 5.43 
6 21.7 rather marked 5.42 
T 23.6 medium 5.42 
8 24.2 very strong 5.41 
9 26.5 strong 5.41 
10 28.1 medium 5.41 
11 30.6 faint 5.41 
12 32.0 faint 5.41 
13 32.4 medium 5.40 
14 34.1 medium 5.39 
15 35.9 medium 5.41 
16 38.6 5.40 
17 39.0 5.40 
18 40.45 5.40 
19 41.55 5.41 
20 44.55 5.40 
21 45.8 5.40 
22 46.8 5.40 
23 48.3 5.40 
24 49.55 51 


Table 3). These rings did not coincide with the additional rings appearing in the pseudo- 
-replicas. | 
The PbS layer evaporated on the film is, in most cases, irrespectively of the kind ` 
of film, composed of crystals (whose sizes vary from 200 to 300 À) which fuse into 
bigger crystals under the influence of illumination by an intense electron beam. 
Photographs 7 and 8 represent a layer not subjected to the action of the intense beam. 
Photographs 9 and 10 represent successive stages of formation of the bigger crystals. 
At first, on the layer composed mostly of crystals whose sizes vary from 200 to 300 À 
there form single, widely dispersed crystals, 1000—1500 À in size. Then in the substra- 
. tum of small crystals, gaps form, the number of big crystals increases and, with a fur- 
ther increase in the intensity of the beam, the big crystals fusg into spheres (Phot. 11). 


Microstructure of Photoconductive Layers 279 


EE‏ ت 


TABLE 8 
No Intensity Plane indices 4525 р 
Ma EY: 
1 faint im 
22 very faint 
3 strong 
4 faint 
5 very strong AT | 6.115 
6 very strong 200 6.175 
7 medium 
8 medium 
9 medium 
10 very strong 220 6.13 
11 faint 
12 strong 31 6.15 
13 medium DEDE 6.12 
14 very faint 
15 very faint 
16 faint 400 6.13 
17 faint 
18 faint 33 1 ο 
19 faint 4 420 6.12 
20 strong 422 6.10 
21 medium 33.3 6.08 
Ë 1 i 
29 faint 440 6.06 
23 medium 
24 medium 442 6.09 
5 
25 medium 620 6.07 
26 medium 622 6.09 
21 medium 551 6.08 
| 1 | 
28 medium 6 04 6.05 
29 medium 642 6.07 
30 nein Ë 5 | 6.07 
(easel 
31 medium ү 4 al 6.07 
820 


Phot. 12 represents another case of changes in the layer due to illumination. On the 
left-hand side the layer has not yet changed and its structure is almost invisible. To 
the right, single crystals having the size of about 300 Å can already be seen. They are 
arranged in chains appearing in the middle of the photograph; finally, on the right- 
hand side of the photograph they form a spongy structure. According to the thickness 
of the layer, the changes due to illumination take different courses. Microscopic 
examination of the layers evaporated دہ‎ the films reveal (though very seldom) places 
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covered with relatively large single crystals having the size of 1000 A whose formation 
is not due to the action of electrons upon the layer, but results from different local 
conditions of condensation, of the layer in the photocell (Phot. 8 and 13). 

In order to make sure that the shift of the diffraction rings in the case of pseudo- 
replicas are due to the change of the structure during the preparation of the specimen, 
and not to improper photocell production, envelopes were made in which plates for 
the pseudo-replicas and specimen holders with the films were placed side by side 
(Fig. 1). These experiments confirmed the results previously obtained, i.e., that at 
direct evaporation on films, a layer giving the proper PbS pattern is obtained, whereas 
the pattern of the pseudo-replica shows characteristic deviations from the structure 
of the NaCl type. 

In the third series of experiments in order to create conditions for which the 
condensation of the layer on the film will more closely resemble the normal conditions 
on the glass, phosphor bronze grids were fixed directly to the glass plate with the 
help of the formvar film covering the grid and the glass surface in its surrounding. 
The results of these experiments did not differ at all from the results obtained for 
layers evaporated on grids placed on the specimen holders. 

The PbS layers evaporated directly from lead tiosulphate (without re-evaporation) 
from the tungsten band have been additionally examined. These layers were homoge- 
neous, within the limits of the resolving power of our microscope. In the diffraction 
patterns there appeared the typical PbS pattern with diffused rings; there were no 
additional rings. 


12 photocells were used for these experiments, from which about 40 specimens 
were examined and 154 photographs made. 


Results 
Diffraction investigations 


1) АП electron diffraction patterns (éxcept those of pseudo-replicas) for layers 
obtained in quite different conditions have a typical PbS pattern. 

2) The PbS structure is classical, regardless of whether the layer is sensitized or 
not (Sosnowski et al. 1947, Piwkowski 1953). 

3) The method of pseudo-replicas seems to change the crystalline lattice of the 
detached layer (either by oxidation with dioxan, or in another way). 

4) Diffraction patterns from pseudo-replicas with regard to the distribution of 
‘intensities and to the ring system suggest, however, that we are not dealing with another 
chemical substance, but with PbS whose lattice was subject to a certain deformation. 
This problem does not seem to be connected with the method of producing the layers, 
but rather with the whole method of preparation. | 

5) The appearance of additional rings should be ascribed to admixtures for cases 
where there is no shifting of the main rings, whilst in the case of pseudo-replicas, 
also to the method of preparation. Also the layers of the photocells gave the PbS 
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Electron micrograph of the 


Phot. 3. 


.  . Phot. 5. Diffraction pattern of a PbS layer, evaporated directly on the collodion film. No additional 
| ۱ À rings are present. 
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Phot. 6. Diffraction pattern of a PbS layer, evaporated directly on the collodion film. 
Additional rings are present 
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Phot. 8. Electron micrograph of a PbS layer not subjected to the action of the intense 
electron beam. Magnification 23,000 X 
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pattern with different (depending on the case) additional rings, or the PbS pattern 
without any additional rings. Therefore, no identification of admixtures causing 
the appearance of these additional rings has been made. 


Microscopic investigations 


1) Microscopic observations allow us to assume that the structure of the layers 
examined under the microscope depend on the thickness of the layers, as well as on 
the conditions of condensation. Here it should be pointed out that even in almost 
identical conditions there appear quite homogeneous domains alongside others covered 
with well-developed, longitudinal or cubic, crystals exceeding 1000 À in size. 

2) It is known from the literature (Cario and Kallweit 1955) that the size of the 
crystals depends on the temperature of the surface on which they condense, on the nature 
of the layer, on the thickness of the evaporated layer, on the kind of substratum 
(physical and chemical properties, cleanliness), on the degree of vacuum, on the pre- 
sence of foreign gases, on the way of obtaining the atom or molecule beam, on the inten- 
sity of the beam and on the time of evaporation. 

Thus there are so many parameters affecting the size and shape of the crystals of 
the evaporated layer, that it is almost impossible to standardize the conditions of the 
experiment so as to obtain a layer of almost identical crystals. 

We may of course determine the general character of the crystals which appear. 
For thin layers we find that the layer is quite homogeneous (the crystals lie beyond the 
resolving power of the microscope), for thicker layers, there appear more often bigger, 
well developed crystals. 

3) The method of pseudo-replicas seems to be indispensable for thicker layers. 

4) In layers obtained in quite different conditions (pseudo-replicas, specimen 
holders with films of different kinds, glass plates covered with formvar film) there 
occur areas where the crystals are big (of the order of 1000 À and more). 

| 5) The microscopic image changes decidedly due to illumination by an intense 
` electron beam. Especially interesting is the formation of crystals from a homogeneous 
layer. 

Changes of PbS layers due to illumination by an electron beam observed under 
a microscope and examined by diffraction will be the subject of a separate paper. 

The authors wish to thank Professor L. Sosnowski for suggesting the subject and 
for valuable discussions concerning the interpretation of the results. 
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Фельтыновский, Глясс, Пивковский, Торунь, Muxpocmpyxmypa фотопроводящих 
слоев сернистого свинца j 


Микроструктура слоев PbS исследована при помощи электронной дифракции и элек- 
тронной микроскопии. Слои получены посредством испарения в вакууме. Препараты для 
исследований сделаны методом формварных псевдо-реплик, а также напыляя непосред- 
ственно на формварные плёнки, коллодийные или алюминиевые. Слои PbS состоят вообще. 
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из кристаллов размера 200—300 А; не замечено влияния рода основания на их величину 
и форму. Дифракционные образы напыленных слоев соответствуют непосредственно 
гранецентрированной решётке типа NaCl. Дифракционные картины псевдо-реплик про- 
являют характерные отступления от этой структуры, заключающиеся B уменьшении диа- 
метра кольца 200 и увеличении диаметра кольца 220 по отношенио к диаметрам остальных 


колец. 
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BEMERKUNGEN ÜBER DIE BERECHNUNG DES POLARISA- 
TIONSGRADES AUS DER KOMPENSATIONSEINSTELLUNG EINES 
GLASPLATTENSATZES 


Von А. ٤ 
I. Physikalisches Institut der Technischen Hochschule, Gdañsk 
( Eingegangen om Januar 1956 ) 


Es wurde gezeigt, dass man auf Grund der Neumannschen Formeln für die Inten- 
sitát des Lichtes nach dem Durchgang durch einen Glasplattensatz (mit Berücksichtigung 
der mehrfachen Reflexionen an den Plattenoberflüchen) eine allgemeine Forrgel für 
den Polarisationsgrad in der Kompensationsmethode von Arago angeben kann. Im 
besonderen erchält man die Formel von Gaviola und Pringsheim, ebenso wie diejenige 
von Sjewtschenko. Auf Grund der allgemeinen theoretischen Formel für eine beliebige 
Anzahl von Glasplatten, kann man also den Polarisationsgrad berechnen ohne mit dem 
Eichen des Satzes zu beginnen. 


E. Gaviola und P. Pringsheim (1924) haben gezeigt, dass bei der Messung des 
Polarisationsgrades eines Lichtstrahls mit Hilfe eines kompensierenden Glasplatten- 
saizes die einfache Formel 

1 — cos** (а — p) 
Р T+ cost (a — p) Q) 


nicht angewandt werden kann. Hier bedeutet k die Zahl der Glasplatten, a und f die 
Einfalls- und Brechungswinkel. Die Formel (1) berücksichtigt nümlich nur den 
direkt durch die Platten hindruchgegangenen Anteil der auffallenden Strahlung, 
vernachlässigt dagegen die mehrfachen Reflexionen zwischen den einzelnen Ober- 
flächen. 

Zwecks Berechnung der Polarisationsgrade die geringer als 0,2 sind, wendet man 
eine kompensierende Glasplatte, dagegen im Falle der Polarisationsgrade von 0,2 
bis 0,5 benutzt man meistenfalls vier Platten. Für eine einzige und für vier Platten 
wurden von Gaviola und Pringsheim (1924) entsprechende Formeln für den Polari- 
sationsgrad angegeben; welche die mehrfachen Reflexionen der auffallenden Strahlung 
zwischen den einzelnen Oberflüchen der Glasplatten berücksichtigen. Auch Sjewtschen- 
ko (s. Lewschin 1951, S. 62.) gab eine regelrechte Formel für den Polarisationsgrad 
für eine einzige Platte. Da es bisher keine strenge theoretische Formel für den Polari- 


(283) 


sationsgrad für einen kompensierenden Glassatz gab (s. Lewschin 1951, S. 62.), musste 
man entweder diesen Satz zunächst eichen oder auf Grund der von Gaviola und Prings- 
heim angegebenen Formeln und Reihen lästige Berechnungen durchführen. 


Jedoch erweist es sich, dass schon im XIX Jahrhundert F. Neumann (s. G. Kirch- 
hoff 1891, S. 166) Formeln für die Intensität des Lichtes, welches von einem System 
von k durchsichtiger Platten durchgelassen und reflektiert wird, abgeleitet hat, 
indem er gerade die mehrfachen Reflexionen zwischen den einzelnen Oberflächen des 
Plattensatzes berücksichtigte. Die Verifikation der Neumannschen Formeln ist von 
Provostaye und Desains (s. E. Verdet 1887, S. 385.) nachgewiesen. Die Übereinstim- 
mung zwischen den Beobachtungen und der Berechnung ist ausserordentlich gut. 


Im folgenden wird gezeigt, dass auf Grund der Neumannschen Formeln man 
eine allgemeine Formel für den Polarisationsgrad des kompensierenden Glasplatten- 
satzes angeben und den Polarisationsgrad für eine beliebige Anzahl von Platten berech- 
nen kann, ohne dass ein Eichen des Satzes nötig ist. Die Neumannschen Formeln 
für das durchgelassene Licht durch einen Satz von k Glasplatten, wenn man mit Zi 
und Г, die Intensitäten des parallel oder senkrecht zur Einfallsebene polarisierten 
Lichtes bezeichnet, lauten 


d D 
lî GE DO ^irgk-n)0-2" (2) 


Ju = 


wobei nach den Fresnelschen Formeln 


r? sin 2a : sin 2p . sin 2a sin 26 
Uem sin? (a + f) cos? (a — В) ' Me sin?(a + В) (9) 


und sin а = n sin B, wo n der Brechungsindex ist. 
Den Polarisationsgrad des Lichtes das teilweise polarisiert ist bezeichnet man, wie 
bekannt, durch 


ml 
ΠΕΠ (4) 


Mit Rücksicht auf (2) folgt dann aus (4) die allgemeine Formel für den Polarisa- 
tionsgrad bei k Glasplatten 


ES k (d — D) 2 
ЖЕЛЕУ - 0804-0 (5) 


Für eine einzige Glasplatte, wenn man іп (5) k — 1 setzt, geht (5) über in 


pain 
ερ гаю (6) 


- was, wie es sich beweisen lässt, mit Gaviolas und Pringsheims Formel übereinstimmt. 
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2 
( E A sin? а 
n 
p= 1 1 =p (7) 
ll 2 [2- 1 + Ce шга | 
n n 


eine mit der von Sjewtschenko identische Formel. 
Für vier Platten (k = 4) unter Berücksichtigung von (3) lässt sich (5) wie folgt schreiben 


2 
٢ ہے‎ a sin? a 
n 


(1 + i |: — (1 + 1. sin? «| — m — sin? o) (n? — sin2 o) 


Wenn man die Formel (8) mit Gaviolas und Pringsheims Formel für vier Platten 
vergleicht, so sieht man, dass die obige Formel viel einfacher ist. Zur Veranschau- 


Neah (3) und (6) bekommt man 


p= (8) 


lichung der Abhängigkeit des Polarisationsgrades vom Winkel a der Glasplattenein- 
stellung für vier Platten (n = 1,52) seien hier die diese Abhängigkeit darstellenden 
Kurven angeführt (Fig. 1.). 


Fig. 1. Polarisationsgrad p als Funktion des Winkels а der Glasplatteneinstellung 
+ berechnet nach der einfachen Formel (1) 
O berechnet nach der Formel von Gaviola und Pringsheim 
A berechnet nach der Formel (8) 
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Die auf Grund von (8) berechneten Polarisationsgrade liegen sehr nahe bzw. sind 
gleich denjenigen, welche man mit Hilfe Gaviolas und Pringsheims Formel erhält. 
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A. Кавский, Заметки на тему вычисления степени поляризации в случае компенсацион- 
ного установления стеклянных пластинок 


Основываясь на формулах Неймана на силу света, проходящего через систему несколь- 
ких прозрачных пластинок, учитывающих многократное отражение от поверхности плас- 
тинок, получено общую формулу, определяющую степень поляризации в компенсацион- 
ном методе для к стеклянных пластинок. В частных случаях получено формулу Гавиоли 
и Принсгейма, а также формулу Севченко. 

Следовательно, при помощи полученной теоретической формулы для любого коли- 
чества стеклянных пластинок, можно вычислить степень поляризации, не прибегая к пред- 
варительной градуировке компенсационного прибора. 
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TWO-CENTRE INTEGRALS IN SOLIDS 


By МАСІЕЈ SUFFCZYNSKI 
Institute of Theoretical Physics, University of Warsaw 
(Received May 7, 1956) 


Fletcher and Wohlfarth's method of calculating the energy two-center integrals 
in solids is described in detail. 


Fletcher and Wohlfarth (1951) using the tight binding approximation calculated 
the energy band of 3d electrons in nickel. In this paper they calculated for the first 
time the necessary energy integrals in the crystal as two-centre integrals between 
the atomic functions centred on the central and the neighbouring atoms. 

The calculation of these integrals will be described here in detail in order to point 
out all the integrals which should be calculated if the Fletcher-Wohlfarth idea were 
to be executed exactly. 

The integrals will be calculated here between the 3d functions only, because these 
have sufficiently pronounced angular dependence to bring out the characteristic 
complexity. 

Furthermore, for 3d electrons in transition metals, the treatment by the tight 
binding method may not be meaningless. 

The integrals will be calculated between functions of o, z, Ó-type symmetry 
with respect to the axis joining the central atom with its neighbour. The tables of 
Slater and Koster (1954) allow the energy integrals in cubic lattices to be found in 
terms of these reduced two-centre integrals: ہے‎ 

The energy integrals in the crystal have the well known form 


Toten ri - U 00 v dex (1) 


The coordinate system originating in the central atom is unprimed, while the 
system originating in the neighbouring atom is primed. 
The 3d functions are 
$ 322-73 
v Gdo) = (те) Fin 
(287) 
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y 4л) = (35) =f, y (847) = | = TOEN 
= 15\3 1513 х2 — у? 
v (3d8) = Ek Zon y (340) = (35) = αρ; 
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The radial part f(r) is taken in the form of a linear combination of Slater functions 
f(D =P? (Ae + Ве"). (3) 


normalized to 1, that is (45/8) (43 a7 + В? b-?) + 1440 AB (a + b) = 1. 
The potential is approximated by (Z is the atomic number) 


V (r)—U (r) = - [1+(Z-1)e”]/ r' for the space outside the atomic sphere of the 
central atom, 


V (r)—U(r) = 0 inside this sphere. 


This is the essence of the idea put forward by Fletcher and Wohlfarth (1951, 1952). 
And it is at this point that the calculation of two-centre energy integrals in solids differs 
from the molecular case; every energy integral must be calculated as the difference. of 
theintegral over all space (оо) less the integral over the atomic sphere of the central 
atom (S). Such an integral will be denoted here without subscript: 


/ = d KS Ji д (4) 
5 

thus (ddo) = (ddo),, — (ddo), etc. for л, ὃ. The radius of the atomic sphere will 

always be taken equal to half the distance between the central and the neighbouring 

atom. Otherwise, the integration over this sphere would result in much too complicated 


formulas. 


Therefore the assumption of Fletcher and Wohlfarth leads us to the calculation 
of the following integrals: 


(ddo) | e 
el = fas (de + Be") E xp Se 6 |ч ew SR ολ х 


(15/4л) (xz) (х =’) 
(15/4л) (ху) (z y’). (5) 


We introduce the spheroidal coordinates as defined by K. Rüdenberg (1951). 


The distance between the centres, that is, between the central atom and the neighbour 
is denoted here by ZR. 


d 16x) (32? —r?) (3272— 2) 
x 


x = x = R[(£—1) (1—72)]* cos p, 0 <p < 2л, 

y = y = R[(&—1) 0—752]^ sin g, | 

z = К(1— ên), lif оо; 

ss 70 — l«mg«l (6) 
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We then have 


r = Rer) 
ο 
& =2А — =, dax = R? (Е? — n?) d£ ہل‎ dg. 


Since in (5) we have to do with linear combinations of Slater functions, we label 
the coefficients 


№ = 

№ = АВ = М, N, = В?, 

MN Ue η, DES ES 

N = AB(Z —1)=N,, (7) 


and the exponents 
On OF == ао 
ar + br' = αι £ + бү} 
ar + (a + ο) r = a, £ + Ban, 
ar + (b + c)r = as Ë + Bam, | (8) 
br + ar = a, Ë + Ва, < 
br + br’ = a; £ + B; n, 
br + (a + c) r = ag + Ben, 
br + (b + ο) г = a, £ + В, 8. 
We.then have 


== d Fe Po = 0, 
a, = R (a + b), ہم‎ — R (b — a), 
а» = R (с + 2a), f he, 
аз = В (а + b + c), ہق‎ = R (b + с — a), etc. 
We can now write 
(ado) = RYN, (ddo) (9) 


and identical linear combinations for the л and 6-integrals. The integrals for a parti- 
cular a and Û are i 


# (dd) 
a متا‎ OU 
(dd) 
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(= 5/8) [9ê — 6(202 ἐπ) FE TF 42272 ο ο 3 
и Е Е ЕСИ 
(— 15/16) [£*n* — 2(£*n* + Em) + £t + η! 4222 7 (0) 


The integrals over all space are well known from molecular studies. The limits of 


1 со 
D = if dn f dé 
оо —1 1 
and the integrals can be expressed in terms of the A, (a) and ,رظ‎ (B) functions of 
Kotani & al. (1940). 


integration are 


(92). = — (8/8) (9 (4,5, AE ο PD 
رھ‎ HAB А,В; АВА, πο... 
CORE + АВ, — АВ, ANE + 9 (А.В, А,В) 
= (SB) [Ba (As = 645 + 94;) + By (—Ay + 64g —949 + В, (—64, +44, ОЬ 
+ B, (64, — 4A,+6A,) + В, (9A; — 645 + 41) + В, (— 9А, + 64. — All, (11) 


aß 
(dam) = (15/4) [By (Ag — A) — В, (A, — 4) — В, (А; — A) F 


+ B; (A, — 49) + B, (A, — 49) — B; (A, — 4), (12) 
ай 
(ddà),, = — (15/16) [(By + B, — 2B,) (4, + A, — 249) — (B, + В, — 280) (A, + 
+ Ay — 245), (13) 
aB 
بت اف‎ ۳ 
3 3 6 8 15 3 3 6 
3 ( او و لے‎ ος, 
15 5 6 15. Te 9 3 
مھ‎ ۹٦۳ (14) 


The last, particularly simple, formula can be used for checking numerical computations. 
The exponent f — 0 gives us the shorter expressions: 


dnt] ta er ope ad dequo rate εἰς 

9 8.2 5,1). نت‎ 
e ЫР (eam 2 ۹ + 75 
×× 3 5 18), ER 


By - (AS ILS ВЕ 
(ao — 85 (14S BUE). (17) 
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The evaluation of the integrals over the atomic sphere constitutes the bulk of the calcu- 
lations. The limits of integration are 


1+7 


ils Јаја 


and the integrals over £ and 7 аге not independent. We use here the abbreviation 


в —a- В. 

ЕВ 5e1 1 4 12 24 120 
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1 2 12 24 168 72 1 23 30 
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In the actual numerical calculations of Fletcher (1951, 1952) and in my own (1955 
a, b) as well, the whole sum (9) was not been evaluated. In the integration over all 
space, only terms with > = 0, 1, 2, 3 in the present notation have been retained, and 
in the integration over the atomic sphere, only terms with у = 0. It was argued that 
а < c < b. For instance, the 3 d electrons in transition metals have atomic functions 
which can be approximated by (3) taking a — 2, b — 5. 

Therefore the formulas for the integrals over the atomic sphere with f #0 
which I have derived (1956) have not yet been used in numerical calculations. 


I thank Dr. G. C. Fletcher for his kind correspondence concerning the two-center 
integrals in solids. 


KPATKOE СОДЕРЖАНИЕ 
M. Суффчинский, Бицентровые интегралы в кристаллах 


Описан метод Fletcher и Waohlfarth вычисления бицентровых интегралов энергии 


в кристаллах. Выписаны формулы для интегралов энергии, вычисленных по атомной 
сфере. 
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SUPPLEMENTARY BOSON-FIELD METHOD AND COLLECTIVE 
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In this paper, the method of second quantization is investigated for an ensemble 
of interacting fermions described by the quantized wave function y. The work here is 
based on the supplementary variables method (D. Zubarev, Candidate's Degree Disser- 
tation, Moscow University, Moscow, 1953). Bohm and Pines, Phys. Rev., 92, 609 (1953). 
А supplementary boson field has been introduced into the theory in such a way that the 
quanta of this field describe the collective motion of an ensemble of fermions (electrons). 
This is possible, after suitable transformations, since ooperatrs connected with the supp- 
lementary boson field appear in the Hamiltonian. It is further shown that the transforma- 
tion of the Hamiltonian can be regarded as a transition to new fermion operators formed 
from the old fermion operators and from boson operators associated with the supplemen- 
tary boson field. The method presented here is compared with the method given in the 
paper of Tomonaga (Progr. theor. Phys. 5, 544 (1950). 


Introduction 


In this paper, the method of second quantization is investigated for an ensemble 
of N mteracting fermions described by the quantized wave function y. 
The Hamiltonian of such an ensemble has the form 
р? 


H= = 1m PH V, y(r) drt 


+ E 1 yt (τι) pt (т) V (In — ral) v (ra) v (г) dr, гу. (1) 


In this formula, m denotes the mass of the fermion, PET V (r), the ο. energy 
of the interaction between the fermions. 

Basing ourselves on the supplementary variables method (Zubarev 1953, Bohm and 
Pines 1953; see also Galasiewicz 1956), we introduce into the theory a supplementary 
“boson field in such а way that the quanta of this field describe the collective motion 
of the group of fermions (electrons). This will be possible, after suitable transforma- 
tions, since operators connected with the supplementary boson field appear in the 
Hamiltonian. There fore we assume that the wave function of the group depends not 


(295) 
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only on the occupation numbers of the fermions, but also on the occupations numbers 
of the bosons associated with the supplementary field. We shall then compare this 
method with the method presented in the paper of Tomonaga (1950). 


Transition to the momentum occupation representation 


In order to put operator (1) in the momentum occupation representation, we 
expand y in terms of the eigenfunctions of the momentum operator: 


y= Dy are la (2) 
1 


Making the assumption that V(r) is a Coulomb potential we then expand the interaction 
energy per unit volume in a Fourier series: 


ἀποὶ ,‏ ,رپ 
V(Ir,—r,) = Y Zoe eik (r= νι) (3)‏ 
k‏ 


(-e is the charge of the electron). In summing over k, we have indicated that the term 
for k = 0 was excluded. This corresponds to regarding the electrons as being on the 
background of the neutralizing ionic charge (see e.g. Bohm and Pines 1952, 1953). 


In the momentum occupation representation, the Hamiltonian of the second quan- 
tization has the form 


hf 1 ‚ 4ле? 
Н = Σ °” α]αι-|- 9 > 2 > aj, aj, Q f, +в Gf, κ. (4) 


Subsidiary condition 


The wave function of the ensemble under consideration depends on the occupation 
numbers of the momentum configuration space and the occupation numbers of boson 
quanta of the supplementary field. Since the boson occupation numbers n, form a set 
of additional dynamical variables, it is necessary to add certain equations. These equa- 
tions will correspond to the subsidiary conditions imposed on the wave function in 
Zubarev's method (1953). In our case, the subsidiary conditions take the form 


(oy + bp = 0, (5) 


where Ûy, bp are the Fourier amplitudes of the supplementary boson field. 


We shall now investigate more closely the restrictions placed on g by imposition 
of subsidiary conditions (5). For this purpose, 


of |... ny... ny: 


we expand the wave function ф in terms 


P= У) Σο а Anka (6) 
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Because of the equality 


ο ipa E= = Vn, + 1 | n, + ln 
bus | Dy, rp, > = үл, | Dy, -k — 1 (7) 


we may write the subsidiary conditions in the form 


21 2 (Vn, +1 «n, п> + Vn, + l <m +l, n 4-1] 9) ny l n ar 


[۱0 ”-k/0 
+ У γην +110, п, ><0,n, +1|9> = 0 (8) 
п. 6/0 
Therefore 
«0,n.,4-1|9 >= 0 
`and 


Mail πρ Vnde Losa danig spl p ο (9) 
for ny, πι κε [0, оо]. 
If we introduce the notation 
<mo ть | o> = C (n. m), ما‎ qu 
then equations (9) can be written 


C (0, m +1) — 0, 


ο. σία ima (11) 
Because of the symmetry with respect to n and m, we obtain, after brief calculation, 
C(n, m) — 0 for n £ т 
and 
С (n. п) = — C (n + 1, n + 1) = C (0,0) (— 1)". (12) 
Hence 
| C (n, n) 2 = | С (0,0) |? (19) 


for п Σο... оо. 


Thus the only states in expansion (6) permitted by the subsidiary conditions are 
those for which the number of bosons with a momentum Ё is equal to the number 
of bosons with a momentum --ᾱ. All these states are equally probable. 


Transformation of the subsidiary conditions 


Our purpose here is to show that by means of the subsidiary condition (5), we can, 
with the help of the unitary transformation 


as سے‎ 
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1 
U = e= exp |- le 7 2 оет (14) 
1 |k|&k, 
reduce Hamiltonian (4) to a form in which appear terms containing the operators 
b, and bjt describing the collective oscillations. The frequency o appearing in equation 
(14) will be determined later when we shall transform the Hamiltonian (4). The 
quantity kç is the quantum mechanical equivalent to the quantity kp = 1/Ap appearing 
in classical plasma theory, where A, is the so-called Debye length (see e.g. Bohm and 
Pines 1952, 1953). 
In this connection we shall now consider the transformation of subsidiary con- 
ditions (5) with the aid of unitary transformation (14). 
The subsidiary conditions can be written 


U (b, + bt,) U*Up = e (b, + b+,) е*Ф = 0 (15) 


where Ф = Ug. 
We shall make use of the following relation which holds for any arbitrary opera- 
tor 4: 


UAU = e Ae = A үү E А] 3i M EYES 


=4+ VY Sa. (16) 


all 


Аз can be seen from (14), S has the form 


ο) ME (DRE او‎ a, sa (17) 


f |к|<һ 
Operators b,, by satisfy the following commutation relations: 
by, by Dp ὂμ = δε 
b, On, — b, b, = O, ' (18) 
by by, — δὲ, ہو‎ 


Employing (16) and taking into account (18), the transformed subsidiary conditions 
have the form 


hk w * گان‎ : 4 
3 Vazh) (0 + 6-а) –е 7 ат-к ارہ‎ 20 (19) 
1 
for |&| < ko. 


We shall make use of these relations when investigating the term in 7 Hamiltonian 
corresponding to the Coulomb interaction. 
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Transformation of the Hamiltonian 


We shall now transform the Hamiltonian (4) using the transformation (14). We 
shall first deal with the part corresponding to kinetic energy. According to (16), 


© ı Го 
ΟΣ я Раза, + Y [^X Λαξ αἰ]. (20) 


If we introduce the notation 


ho 


» 
| τ e = (БЕЛЛ, e (21) 
(а, commutes with the operators αχ), then 


SER ууа аа (22) 


kISk f 
Operators a, satisfy the following commutation relations: 
+ a+ € 
Gy, ак, + Gk, Gk, = On, ko 
αν, αν, + ak, % = 0, 
de di |- ας ay = 0, 7 2 


Making use of the above, we obtain 


Ë Y fat οι] = x 3 а, k 5 SS d а} ка, 
f 


f |kl<k 
(2) 
Ё 2 Zë аў 7 =2 M 9 ок, αι (R) ko) а k-k, ау, 

1 Í EE 

E Goen ap | 0 (24) 
" αι αι == 39 
f 
Therefore (n) 
Ё NATI «j| = 0 for n 3 (25) 
1 


Thus the series of terms in the sum in equation (2) terminates after the second term. 
We shall now take up the Coulomb interaction. This can be split into short-range 
and long-range parts (k < ko, k > ko): 


1 4ле? TL T4 dh 
9 5. 7 Qf, Qf, ау, +в Gf, к 
4 k ` ї fa 


| 4ле? 
=Ë Σ tg Ў | ey οἱ Gf, й+к аук. (26) 


|К| < ko |k| >ko f» В 
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We shall represent the Coulomb interaction in another form. This will be done 
for two reasons: firstly, in order to utilize subsidiary condition (19) in a simple manner, 
and secondly, in order to show that the term corresponding to the Coulomb inter- 
action does not undergo any change in consequence of the application of transformation 
U to the Hamiltonian (4). 


If we make use of the relations 


T E uu 
OF ور‎ Ope Opn = ар, αι κα арав Gf Of tbe (27) 
and 
ات‎ MR 
Ха а= № 
1 
then the long-range interaction can be represented in the form 
, E + 1 ‚ 4ле? 
ΣΣ - > aj «x f, 3 ака – 5 У" پر‎ (28) 
(ap? f. k <ko 


to which subsidiary condition (19) is easily applied. Since we have the following equality 
holding for the operators (operating on wave function ®): 


X hk | о \* ü 
> a ر‎ = -z (aog) (be + b3,) for | k | < ko (29) 
we can then replace the sums >) af ,, a; and >} رہ‎ και in the first term of equa- 


1 


tion (28) by the corresponding expressions containing the boson amplitudes b, and ` 
bi. This term then becomes 


)30( 7720ھ 


[К] >4 


From the right-side of equation (24.2), we now separate out the terms of which 


= = i> (be — bk) (bx — b_k). (31) 


[| < № 


4ле? N 


If we assume that ὦ = (frequency of the Langmuir oscillations of the 


` plasma), then (31) reduces to the qud 


о. (32) 
|к|<к, ү 


and сап be joined to (30). Аз a result, we obtain part of the transformed Hamiltonian, 
which we denote by Н 


osc* 


= Ὁ ho (bi b, + и) (33) 
|| Κο 
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This is the Hamiltonian of the second quantization of oscillators vibrating with the 
Langmuir frequency. It describes the collective oscillations of an ensemble of fermions. 
(It can be shown that assuming a Coulomb interaction between fermions, the separation 
of Н is possible only if w is the frequency of the Langmuir oscillations.) 

The part of the Hamiltonian corresponding to the Coulomb energy commutes 
with transformation U. This can be most easily shown by bringing both terms of (26) 
to the form given in equation (28). Thus, after employing the subsidiary conditions, 
Hamiltonian (4) transformed with the help of (14) takes the form 


р? f? 4ле? 
H'— UHU* = Y “аў ay + У ho (bit δι + V) V У N+ 
z μάζα, kh, 
Zahl? (hk 
ολο ыы s 
kish f 2 
е? 2zh + ا‎ 
ys 2m X Σ ο (Er, в) (b-m — bi.) (bi, — br,) ور ور‎ αι + 

ër kl<k f 

Κο| < Ro 

kı +k, +0 


1 4ле? 
T 9 » XY P af, Of, Of, µε а-к 


|k| > Κο Її, ГА 


where ἐμ is a unit vector parallel to vector А. This Hamiltonian describes interacting 
fermions by means of a short range (screened — see Bohm and Pines 1953) Coulomb 
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potential, bosons with an energy hw | ὦ = and an interaction between 


the bosons and fermions. (In Zubarev's paper (1953), Hamiltonian (34) was obtained 
by a method different from that given here.) 


Transition to new fermion operators 


The transformation of Hamiltonian (4) with the help of (14) can be regarded as 
a transition from operators a, to new operators پا‎ Thus Н’ = UHU* can also be 
written in the form 


AA 1 ' Але? 
H'= Oe E Af Arts 5 FEY № Aj, Aj, Anin Аһ (35) 
1 k fo fs 
where 4, = Ua, U+. 
Using equation (16), we obtain 


Е 


55 n un 
A, = Ua,U* = ау = [5 а, (36) 
n 


ο... ;‏ ہے 
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Taking into account (22) and (23), we calculate the respective commutator brackets 
(1) 


[S, αι] = > Opn аук, " (37) 
|k|< ko 
(2) 
Sa, = а, Ak, αι. 
[ d E kı ^k. “fH ki + Κε Я (38) 
к. < Ро 


Ву induction, it сап be shown that 


(п) 


[S, αι] = ὍΣ ак, ἄμ, ... Ok, Oft ky ped + ker (39) 
К | S ko, 


[Ку | S ko 


Hence, after replacing а, from (21) we obtain 


со ΤΉΝΕ Sup) 
^ (b-k, — бк,) (b-r, κ, 
AT а; + уу ОЕ ... pp ad и’. D 
ni <А μοι, : 2 A 
See? (δι, -- bi) (bn Б), 
Af = αἱ + Y ( 24 „4 ο, چوس ہر ریت‎ а; ўњ; (40) 
ai dh “|k,|<k 1 " ш 
where f — A e is a small parameter. 
w 


U is a unitary transformation, and therefore, operators A,, like operators ау, 
describe a fermion field. From the point of view of the operator (35) acting on 
the wave function Ф = Uy operators A,, А} are „old“ and operators a, دو‎ „new“ 


Comparison with Tomonaga's method 


The use of Zubarev's and Bohm and Pines idea of supplementary variables 
(1953) for investigating a group of interacting fermions on the basis of a second 
quantization gives us a new possibility for comparing the elementary perturbation 
method of the plasma oscillation type with the method proposed by Tomonaga (1950). 

Since the operators 4 ү describe a fermion field, the “fermion branch“ in the energy 
spectrum is conserved. The new operators are constructed from fermion operators ap 
and from boson operators b, associated with the boson field. As a result, there appears 
in the Hamiltonian a quanta energy corresponding to the collective oscillations. This 
problem was solved for the three-dimensional case. 

In the method proposed by Tomonaga, a transition is made from operators a, 
which describe a fermion field, to operators ор, which describe a boson field. Here 
the “fermion branch“ does not appear in the energy spectrum. The boson operators op 
are constructed from the fermion operators а». The problem was solved for the one- 
-dimensional case. 


I would like to thank Professor M. Günther for helpful discussion. 
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З. Галясевич, Метод „добавочного“ бозонового поля, а коллективное колебание 

В работе изучается посредством метода вторичного квантоваыия система воздейству- 
ющих фермионов, описанный сквантованной волновой функцией у. По примеру метода пере- 
„добавочных“ менных (Зубарев Д. Кандидатская дисеертация, Московский Университет, 
Москва 1953 Бохм и Пинес, Phys. Rev., 92, 609 (1653), введено в теорию „добавочное“ 
бозоновое поле так, чтобы кванты этого поля описывали коллективное движение системы 
фермионов (электронов). Это возможно после соответствующих преобразований, когда в 
гамильтониане появятся операторы, связанные с „добавочным“ бозоновым полем. Далыпе 
доказано, что преобразование гамильтониана можно считать переходом к новым фермио- 
новым операторам, построенным из старых фермионовых операторов и операторов 60- 
зоновых, связанных с „добавочным“ бозоновым полем. Представленный метод сравнен 
с методом в труде Tomonaro [Progr. Theor. Phys. 5, 544 1950)]. 
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This paper discusses a possible explanation for the diurnal changes in the intensity 
of the green line of the night glow. These changes may be evoked by periodic changes of 
pressure in the upper layers of the atmosphere arising from the tidal activity of the sun and 
the moon.The correctness of the above hypothesis can be confirmed by investigating 
the relation between the position of the moment of maximum green line intensity and 
the lunar phase. It can be seen from the curve drawn on the basis of published curves of 
green line intensity changes available to the author that such a relationship actually exists. 


Photometry of the night glow spectrum provides very interesting results. Among 
those facts which have been established is the occurrence of diurnal variations in the 


intensity of the green line 5577 À arising from the forbidden transition 15, > 1D, 
of neutral oxygen. The character of these changes is rather curious and so far has not 
been fully explained. The intensity of the green line increases during the night, generally 
achieving a sharp maximum in the interval between 23 A. and 03 ۸۰ local time, and then 
decreases. Up to the present, there has been no explanation in the literature on this que- 
stion of the changes of the position of the maximum green line intensity during the 
night. Some possibility of explaining this phenomenon is provided by considering the 
general conditions under which the forbidden transitions occur. The green line arises, 
as mentioned above, as a result of the transition of the atom from the 15, state to the 
1D, state. They are both metastable states of the oxygen atom. The oxygen atom, 
excited to the state 15, can remain in this state a relatively long time (mean life is 0.5 
sec) and make the transition to one of the lower levels. The atom may leave the meta- 
stable state by: 

1) dropping to a lower energy level with the radiation of a quanta, ` 

2) absorption of radiation — the atom makes a transition to a higher energy level, 

3) a collision of the second kind with another atom or molecule, 

4) a collision with an electron — the excitation energy of the atom is entirely 
converted into the kinetic energy of the electron. 


(305) 
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The probability of the transition from a metastable state by means of a collision 
of the second kind with atoms and particles is, in general, very great when the gas 
density is high. This explains the fact that it is very difficult to obtain the forbidden 
radiation in pure gases under laboratory conditions. The excitation energy of the atoms 
or molecules causing the transition from the metastable states must be less than 
the energy for exciting the atom to the metastable state. If the excitation energy 
of the lowest excited level of the atoms causing the transition from the metastable state 
is greater than the energy of the metastable state, an elastic collision occurs without 
the loss of energy. For example, the oxygen atom can remain in the metastable state 
in an atmosphere of inert gases, since their lowest excitation energy is several times 
greater than the energy of the metastable state. Use is made of this in obtaining forbid- ` 
den radiation under laboratory conditions. If there are present atoms which can cause 
transitions in atoms in a metastable state, radiation of the forbidden lines can occur 
only when 


τε <? (1) 


where т is the mean time between two such collisions, 9 is the mean lifetime of the meta- 
stable state, and & is the efficiency coefficient for the collision. The mean number 
of collisions which suffices to cause the transition from the metastable state is given 
1 es $ v : 
by the value of — . This is a first approximation of the problem; it assumes that all 
€ 
atoms have the same lifetime. In reality, many atoms have a shorter lifetime in the meta- 
stable state. i 
Transitions from metastable states of atmospheric oxygen are probably caused 
by molecular nitrogen and oxygen. These molecules have densely distributed energy. 
levels very near the metastable energy levels of oxygen. In this case can be taken roughly 
as equal to unity. 
It is well-known from molecular physics that 
T= x -—— I (2) 
р 
where p is the gas pressure, T'is the temperature, and x is a constant for the given gas. 


At constant temperature, т depends only on p. Аз а result of condition (1) the forbidden 
radiation can occur only from a certain value: 


үт 


pa < x ہے‎ (3) 
If we consider the possibility of forbidden radiation arising in the earth's atmosphe- 
re, we see that this can occur in those layers beginning at a height h, corresponding 
to a pressure p,. 


For the sake of simplicity, we assume here that the temperature of these layers 
is constant and equal to 240°K. 
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We shall calculate the intensity of forbidden radiation falling within a unit solid 
angle. Small zenith distances will be taken for simplicity. It can be seen from Fig. 1 
that the intensity dI arising from an element of volume at а height A is proportional to 
the number of excited atoms in unit volume n (A), to the size of the element and inver- 


Fig. 1 


sely proportional to the square of the distance from the observer O. The element of 
volume is cut out from a layer of thickness dh by a unit solid angle and has a base 
proportional to A? and height dh. We can write, 


h) h? dh i 
dI = ЖУ ч =an(h)dh b. (4) 
where a is a proportionality factor. 
The complete intensity will be 
1 — a | n(h) dh (5) 
hd 


For heights greater than Р, the number of excited atoms in a unit volume will be pro- 
portional to the total number of oxygen atoms under consideration. It can be assumed 
that the number of excited atoms changes in accordance with the function 


r. (b) = مد‎ fb (6) 
Substituting this into the expression for [ d 
со 
I (ha) = ang ТА e fh dh = وہ‎ e Pha (7) 
hg ñ ү! 


We thus see that changes in intensity of the forbidden lines will take place when 
changes of pressure, and therefore of height ἦα, occur in the upper layers of the 


atmosphere. 
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It so happens that such changes of pressure in the earth's atmosphere actually 
occur. From the theory of atmospheric tides it is known that as a result of the gravi- 
tational effects of the sun and moon, tides arise in the atmosphere. Maximum rare- 
faction occurs at the time of the lower culmination ofthe sun and moon. During the 
new moon the effect is additive. On the surface of the earth this tide phenomenon is 
very weak and therefore it is difficult to establish. The amplitude of the solar tidal wave 
is 0.26% of the barometric pressure. The activity due to the moon is still weaker. 
The amplitude of the lunar wave is barely 1/16 of the solar wave amplitude. The pheno- 
menon is quite distinct in upper layers of the atmosphere. It is known from Appleton 
and Weekes’ investigation of the behaviour of radio waves in the E region of the iono- 
sphere that at this height (above 100 km), the amplitude of the lunar wave is about 
10% of the static pressure. If we denote the change of pressure resulting from the tide 
(caused by the moon) by ما‎ and the static pressure by po, then the relative change 


of pressure will be 
(22) = 104 ES (8) 
Ро / h= 100 km Po | h-0 


There the solar amplitude is greater by nearly one order of magnitude. | 
Let us return to the calculation of the effect of a change in pressure on the intensity 
of the forbidden line. For the increase in intensity we obtain the equation 


4 I (ha) = αἱ Aha = — аще Aha (9) 
dha 
The relative change in intensity will be 
I ۱ 
2 = — BAhg (10) 


The quantity Ah; can be calculated if we know the change of pressure Др at height 
ha, since ! 

p (1) = poe” ` (11) 
(for the earth's atmosphere y = 1.5 x 10-4 em!) then 


1 
Ah; = yp (ha) Ap (ha) (12) 


These formulae connect the relative increase in intensity of the forbidden line with 
changes of pressure in the upper atmosphere. 

We see that in considering the conditions giving rise to forbidden radiation from 
oxygen in the atmosphere, the density changes resulting from the tides cannot be ig- 
nored. Around the time of the new moon the tidal waves add to each other; maximum 
rarefaction, and therefore, maximum green line intensity, does not however occur 


at midnight, local time, but about one hour afterwards (analogous to the port time in . 


marine tides) as a result of the inertia of the atmosphere. After the new moon, the maxi- 
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mum of the lunar wave occurs later, and before the new moon, it occurs earlier. The 
resultant maximum rarefactions and the maximum green line intensities which 
accompany them should be similarly distributed in time. The photometry of the night 
sky spectrum can be made only around the time of the new moon, since, in this period, 
the light of the moon will cause no interference. The correctness of the above expla- 
nation can be checked by investigating the correlation between the position of the maxi- 
mum green line intensity and the lunar phase. With this aim in mind, the times of 


MP (h) 


BOE CE οἱ ο. μα d 
о - Elvey & Roach Mn s (d) 
е – Barbier, Dufay & Williams 


Fig. 2 


maximum green line intensity for different nights were taken from curves published 
by Elvey and Roach, and Barbier, Dufay and Williams; the corresponding times of 
the new moon were taken from issues of the „Nautical Almanac“. Denoting the time 
of the maximum by M,, the time of the new moon by M, local midnight by P, 
the time interval М, — P,, in hours and M, — P,, in days, were calculated for the 26 
curves that were used. I used all the curves published by Elvey. Of those published 
by Barbier (where the accuracy of the measurements and of the reading of the maxi- 
mum point was considerably less) only the curves with a distinct maximum were used, 
Some curves (especially Elvey's) had a secondary maximum weaker than the main 
one. Їп these cases also, the times corresponding to the main maximum were taken. 
The data gathered cover the period running from the fourth night immediately pre- 
ceding the new moon to the fifth night immediately following .The times of the maxima, 
and therefore th: intervals M,,— P,,are determined to an error of + 0.5 h. The curve 
shown represents the time interval M,, — P,, as the function of the lunar phase, and 
therefore of the time interval М, — P,. The appearance of a distinct dependence 
between the moment when maximum green line intensity occurs and the lunar phase 
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seems to support the correctness of the suggested explanation. The discrepancies which 
appear can be explained not only by errors of measurement, but also by the dependency 
of the magnitude of the tides on the distance between the earth and the moon at the 
time of the new moon; this distance is known to be variable. For a final confirmation 
of the hypothesis presented, and of the curves obtained (Fig. 2), it is necessary to make 
a considerably greater number of observations. 

In conclusion, I would like to thank Professor Н. Niewodniczanski for his valuable 
remarks on the problem of forbidden radiation and for his helpful discussion on the 
above problem. 
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JI. Лишка, О некоторой возможности объяснения суточных изменений напряжения 
зелёной линии свечения ночного неба ОТ 5577 А. 


Настоящая работа занимается некоторой возможностью объяснения суточных изме- 
нений напряжения зелёной линии свечения ночного неба. Они могут быть вызваны перио- 
дическими изченениями в высоких слоях атмосферы. Происходят они от приливного дей- 
ствия солнца и луны. Справедливость вышеупомянутой гипотезы можно удостоверить, 
изучая зависимость между положением момента максималной интенсивности зелёной линии 
и фазой луны. Из диаграммы, совершенной на основании доступных автору опублико- 


ванных кривых изменений интенсивности зелёной линии, видно, что такая зависимость 
существует.в самом деле. 
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HALBKLASSISCHES BILD DER SPINWELLE IM FERRIMAGNE- 
TIKUM 


VON H. CorTA 


Laboratorium für Ferromagnetika des Physikalischen Instituts der Polnischen Akademie der Wissenschaft, 
Poznan 


( Eingegangen am 4 Juni 1956) 


Im ersten Teil ist eine kurze Begriindung des halbklassischen Spinwellenbildes 
angegeben. Dabei wird eine Differenzengleichung fiir die lineare Atomkette gebraucht, 
anstatt der oft benutzten Differentialgleichung fiir die kontinuierliche Spinverteilung. 

Im Hauptteil sind die durch Keffer, Kaplan und Yafet auf Grund des halbklassischen 
Spinwellenbildes für Ferromagnetika und Antiferromagnetika hergeleiteten Dispersion- 
formeln о (k), auf den Fall des Ferrimagnetismus und Antiferrimagnetismus im Sinne 
von Smart verallgemeinert. In der Smartschen Nomenklatur bedeutet der Ferrimagnetis- 
mus gleichgerichtete Nachbarspins, und der Antiferrimagnetismus — entgegengesetzt 
gerichtete. In beiden Fällen ergibt sich eine gleichartige Abhängigkeit o (k) نہ‎ k?, die 
mit dem makroskopischen Verhalten beider Magnetikumstypen übereinstimmt. Die 
Formel w (k) ~ k für Antiferromagnetismus lässt sich nicht durch gewöhnlichen Grenz- 
übergang S, > δι, d. h. durch Betrachtung des Antiferromagnetismus als speziellen Fall 
des Antiferrimagnetismus, erhalten. 


I. Die Spin-Präzessionswelle 


Spinwellen treten bekanntlich schon in der Theorie auf, die auf dem Heisenberg- 
schen Modell gestützt ist. Die physikalische Deutung der Lósungen des Bloch-Slater- 


schen Gleichungssystems führt nämlich zu einer Spinumkehrungs-Welle. Ein genaue- 


res Bild dieser Erscheinung erhält man, wenn man aus der quantenmechanischen Ве- 
wegungsgleichung für den m-ten Spin ausgeht (die sogenannte phänomenologische 


Auffassung, s. Heller u. Kramers 1934) ο; 


ih Eom [Ss gal, = (1) 


wobei % der Energie-Operator des Systems und > der Spindrehimpulsoperator für 


das m-te Atom ist (hier hat er die Dimension des Drehimpulses). Wird die Austausch- 


energie durch die Diracsche Formel (Dirac 1947) ausgedrückt, dann ist die Gesamt- 
(511) 


| 
$ 
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energie % in einem konstanten magnetischen Feld H, gleich 


й=-у E Ba уу, 45 (2) 


wobei у = e/m der giromagnetische Faktor für Spinmomente und 4;; das Austauschin- 
tegral zwischen dem ;-ten und j-ten Spin bedeutet. Setzen wir Ф in die Gl. (1) ein, 
so lässt sich der Kommutator zur rechten Seite auf Grund der Vertauschungsrelationen 
umformen (siehe Herring и. Kittel 1951, Appendix A) und dann ergibt sich 


4 > 2 5 
یں‎ و٦‎ H 
J 


Der Übergang von den Operatoren zu den Erwartungswerten ist mit dem Über- 
gang zum halbklassischen Bild gleichbedeutend. Dieses Bild ist selbstverständlich 
vom streng quantentheoretischen verschieden. Z. B. im Falle der linearen Atomkette 
haben wir in halbklassischer Näherung die Gleichung (wenn wir nur Austauschwir- 
kungen zwischen den Nachbaratomen berücksichtigen) 


dŠ, ë = 2A > ج‎ 
i = Kë X μον H, + FX] oc S 5. ο] (4) 
Ihre physikalische Deutung ist leicht ersichtlich, wenn wir zur rechten Seite ein Null- 
glied | 
ας | 
м9, Ж Ἔ S5 (5) 
hinzufügen und erhalten so 
dore وت‎ EE E 
ΠΕ = Sin X d ^ H, ar he IG 59 nn! s (6) 


Diese Gleichung ist vom Typus einer Präzessionswellengleichung, weil das Glied 
Sm+ı — $9 = (Su 51) 
eine zweite Differenz darstellt und damit der zweiten Ableitung entspricht. Die strenge 


quantenmechanische Gleichung ist nicht von demselben Charakter, da fiir die Opera- 
toren der Ausdruck (5) infolge der Unvertauschbarkeit von Spinkomponenten nicht 
verschwindet. > 

Die exakte Lösung von Gl. (4) bzw. der äquivalenten Cl. (6) ergibt eine Spin- 
-Prüzessionwelle 
Sm = Sms + В cos (t) + C sin фи (0) OH 
mit g, (t) = akm — ot 
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Dabei muss 5, | B | C und | B Ces C| gelten. Die Konstante a bedeutet hier den 
Abstand ' Nachbaratomen. Die Wellenzahl E ist mit w mittels der Beziehung 


w = us y Hy Ree 75 N — cos ak) Bas (8) 


gebunden.Der Ausdruck (7) ist eine Lósung der Gl.(6) nur dann, wenn S, , die vektorielle 


Komponente von Se um die durch H, festgelegte z-Achse ist. Dieser το H, setzt 
also die Richtung der Präzessionsachse fest und die Winkelgeschwindigkeit der 
Präzession können wir wie folgt schreiben 


— 


ЛА = 
о = + μον Hy + za (1 — cos ak) Sal : (9) 


Die Richtung der Präzessionsachse ist für sämtliche Atome dieselbe. Die Rolle des 


konstanten Vektors H, wird von der Resultante des äusseren Magnetfeldes und des 
Anisotropiefeldes gespielt. Auch ein drittes Glied kann mit inbegriffen werden, näm- 


lich ein Ersatzfeld H, für den Einfluss von Leitungselektronen. Nach Heber (1952a, 
1952b, 1953) lässt sich dieser Einfluss nüherungsweise dadurch beschreiben, dass 


in den Klammern auf der rechten Seite von (3) noch ein weiteres Glied x S, erscheint, 
das zum Gesamtspin S, aller n Atome proportional ist. Der Zustand mit einer Spinwelle 


ist der absoluten Sättigung so nahe, das der Vektor Sọ als konstant betrachtet werden 
kann, und somit 


dp об 


E x > 
H, 


Da H, als Summand des Vektors H, auftreten kann, dürfte sich die Beteiligung уоп, 
Leitungselektronen in dem Betrag des Larmorschen Gliedes der Winkelgeschwindig- 
keit äussern, natürlich in sehr tiefen Temperaturen. Die Präzessionswelle schreitet 
also nach (7) in solcher Weise fort, dass die Phasendifferenz O zwischen Nachbaratomen 
gleich 


9 = νι — Фи = ak (10) 
e : 
ist. Für kleine Wellenzahlen £ < (1/а) wird 1 — cos ak = 5 a? k2 und (8) geht in 
die bekannte Formel 
24 = 
و‎ = ру Ho + ++ а? K Smz (11) 


über, die sich aus der Präzessionswellengleichung für das kontinuierliche Spinvektor- 


feld ergibt (siehe z. B. Döring 1948). 
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II. Spinwellen in einem ferromagnetischen und antiferromagnetischen Medium 


nach Auffassung von Keffer, Kaplan und Yafet 


Zusammenfassend kónnen wir sagen, dass das halbklassische Bild der Spinwelle 
sich auf dem Vektormodell eines um eine feste Achse präzessierenden Spins stützt. 
Im Grundzustand ist für jeden Gitterpunkt [οἱ = Baas worin Spaz den maximalen 
Betrag der سو ا‎ рова Für den ersten angeregten Zustand 
gilt JS = — (1/п) (siehe Heller и. Kramers 1934). In den hóheren Zustünden 
wird E =. к vermindert, weil jede weitere Spinwelle die |5„„| um Lin verklei- 


nert; der Prüzessionsradius 
R = VS: +5 


nimmt dadurch mit der Temperatur zul. 

Unter Anwendung dieses einfachen geometrischen Modells, haben Keffer, Kaplan 
und Yafet (1953) auf elementarem Wege die angenäherte Abhängigkeit o (k) für 
Ferromagnetika abgeleitet?, nämlich 


w = ug y H, — Gars (11a) 


Der Unterschied mit (11) besteht nur darin, dass anstatt der Komponente S, der abso- 
lute Betrag,S des Spinvektors auftritt. (Das Vorzeichen des Austauschgliedes stimmt 


mit (11) überein, weil S gegenüber H, die entgegengesetzte Richtung hat. Auf 
ühnliche Weise haben die obengenannten Verfasser auch die analoge COS für 
Antiferromagnetika erhalten, nämlich 


ὦ ~ iy y شر‎ ak (12) 
die schon früher auf anderem Wege abgeleitet worden ist (s. Hulthén 1936 und Ander- 
son 1952). Dabei hat es sich erwiesen, dass bei abwechselnd gerichteten Nachbarspins 
eine Präzessionswelle im strengerem Sinne (d. h. w wird für alle Atome gleich) nur dann 
durchlaufen kann, wenn der Präzessionsradius R, für Rechtsspins von dem Präzes- 
sionsradius R, für Linksspins verschieden wird. Infolgedessen werden auch die Kompo- 


nenten S, für ΤῊΝ beiden Spinrichtungen verschieden und damit wird eine resultierende 
Komponente Xs, „ 0 des ganzen Gitters auftreten. Das Zeichen dieser Resultante, 


wie das Keffer (1952) erwiesen hat, hängt nochvom Wellentypus ab (es kann nämlich 
entweder R, < R, oder R, > R, sein), und ihr Betrag nimmt für jede erregte Spin- 
eo E 


1 В hat die Dimension eines Drehimpulses. 


? Dazu genügt eine Betrachtung der linearen Kette. Im See “ΠΝ wir uns nur auf diesen 
Fall beschrünken. 
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welle desselben Typus um Ё zu. Die Gleichberechtigung beider Wellentypen führt 
zu Antiferromagnetismus. Bei Überlegenheit nur eines Wellentypus haben wir dagegen 
Ferrimagnetismus, obwohl alle Spinvektoren gleichen absoluten Betrag haben. 


III. Spinwellen in einem beliebigen Ferrimagnetikum 


Die Betrachtungen von Keffer, Kaplan u. Yafet kónnen für den Fall verschiedener 
Beträge S, = S, benachbarter Spinvektoren verallgemeinert werden. Es wird bequem 
sein hier die Nomenklatur, die Smart (1955, S. 369) vorgeschlagen hat, zu übernehmen. 
Mit dem Namen „Ferrimagnetismus“ werden wir nämlich den Fall bezeichnen, wo 
die benachbarten Spinvektoren gleichgerichtet sind, mit „Antiferrimagnetismus“ 
dagegen einen solchen mit entgegengesetzten Nachbarspinvektoren. Noch allgemeiner: 
wir werden über „Ferrimagnetismus“ reden, wenn die Spinvektoren in verschiedenen 
Gitter-Überstrukturen gleichgerichtet sind, über „Antiferrimagnetismus“ dagegen, 
wenn diese entgegengesetzt sind. Mithin werden wir den Ferrimagnetismus im Néel- 
schen Sinne als Antiferrimagnetismus bezeichnen. Bei einer solchen Nomenklatur 
erweist sich der Ferromagnetismus als spezieller Fall des Ferrimagnetismus für 
Sı = S, und der Antiferromagnetismus als spezieller Fall von Antiferrimagnetismus 


für Ki = Ka 


a) Der ferrimagnetische Fall 


Auf ähnliche Weise wie bei Keffer, Kaplan und Yafet und mit denselben Nähe- 
rungen werden wir die Beziehung @ (k) für das ferrimagnetische Medium im Sinne 
von Smart herleiten. Einfachheitshalber setzen wir H, = 0, da das Larmorsche Glied 
in allen behandelten Fallen identisch bleibt. 


Bezeichnen wir die Summe benachbarter Spinvektoren Sp ME S41 durch δη; wenn 


_ beide Vektoren zur Überstruktur I gehóren, und durch و‎ wenn sie zur Überstruktur 
II gehören, dann ist nach Gl. (4) das auf den m-ten Spinvektor wirkende Drehmoment 


für die Überstruktur I 


Οι = 21 δι 2 sin Ф, (13а) 
und für die Überstruktur П 
G, = — E وک‎ Zu sin $, (13b) 


gleich (S. Abb. 1.). Wegen der Kleinheit von 0, Фу und Ф, können wir 
Σι-- 25, | | 242 == 25. 


х А "d 7 
sin Фу == $ E == 5, == 5 (14) 
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setzen. Für die Prüzessionsbewegung gilt 


G 
w = R (15) 
daher führt die Forderung 
оу = 0 = و‎ | (16) 
zu der Bedingung 
Mit Hilfe von elementaren geometrischen Überlegungen berechnet man (siehe Keffer’ 3 
Kaplan u. Yafet 1953, Fig. 1) 


EC H, = ОМ 


s = 0с SÉ RT 
Sh mx OF ہے‎ I GF "3 wa 
ig^ ابس پوت‎ р qo er sai مم زط طرف سنہ‎ 
فو‎ deba رام سم تقر‎ ac ades جازم‎ iom οσα 
ae jd MI rr 
gy. WI BOE „с Cie E Li NES 1 айкай 095115670 nay 


MER 
μις et! Aas nU gg me. amd ond E eeng ` zeit? yaz 00 
πο Geer " 
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und schliesslich 

A 44 3i: وک‎ 
h? S, + 5, 

Für S, = 5, = S geht die abgeleitete Formel in die Abhängigkeit (11a) mit Н, = 0 über 


а? k? (22) 


b) Der antiferrimagnetische Fall 


Die analoge Beziehung für das antiferrimagnetische Medium wird im Folgenden 
an Hand von Abbildung 2 hergeleitet. Mit den oben angeführten Bezeichnungen 
ergibt sich 


2А ч 
G, = — Τα δι У, sin (л — Ф,) für Uberstruktur I, 
24 А жү 
С. = P S, Ул sin (z — Ф,) für Uberstruktur II (23) 
und mit den Näherungen (14) 
4A 
4A ~ 
Es gilt nun ۱ 
کے‎ —— _ ба 
x gon C xis s (25) 
mithin ergibt die Bedingung (16) 
Sr Ж 26 
SE (26) 


Ahnlich wie im vorhergehenden Falle kann man 


ак 2, 


x = Rp Rp ee; 
; a? 12 
рыл Se рЫ ہے‎ I (27) 
erhalten (siehe Keffer, Kaplan и. Yafet 1953, Fig. 2) und daher 
_ а? k2 R, R + Ri 
К. — R, SC NIE (28) 
Für kleine Ё wird also 
ος EM edere ceri 
_ atk? R, ور‎ + Ri (29) 


AR Be 
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Mit den Näherungen (20) erhält man 


πῶς 
رھ‎ s ΑἹ + 1 Тоя Sy 


۰ = (30) 
2 Og š 2 р? 4 
erai 

S) ni- αἴ 


K" 


Abb. 2 Achsenschnitt durch die betrachtete Prazessionskegeln im Falle des Antiferrimagnetismus 


706 р “ЗС H, — MO 
5) = OF’ ОЕ! = FO Ro CP 
TED R = EA 
pH NER CE 
У, OK = OK" — ОК y = DA 
und schliesslich 
سنا مر‎ 
رہ‎ = ————*.gi[ 31 


Um die Formel (12) zu erhalten, muss man schon in (27) К’, = R, setzen und dann 
Keffer, Kaplan u. Yafet folgen. 
Herrn Prof. Dr. Szczepan Szczeniowski fühle ich mich verpflichtet, für die mit 


mir durchgeführte Besprechung der obigen Arbeit und die genaue Durchsicht des 
Manuskriptes, meinen besten Dank auszusprechen. 


KPATKOE COJIEPXXAHHE 


Г. Цофта, Полуклассическая картина спиновой волны в ферримагнетизме 


В первой части настоящей работы дана в сокращённом виде мотивировка полукла- 
ссического изображения спиновой волны. Использовано при этом дифференционные уравне- 


для для линейной цепи спинов вместо применяемого часто дифференционьного уравнения 
ние непреривного размещения спинов. 
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В коренной же части работы дисперсионные формулы o (k), полученные, на основании 
полуклассического изображения спиновых волн, Кефером, Капланоми Яфетом для 
ферромагнетиков и для антиферромагнетиков, обобщены на случай ферримагнетизма 
и антиферримагнетизма в толковании Смарта. Ферримагнетизм в номенклатуре 
Смарта обозначает согласованно направленные соседние спины, а антиферримагнетизм — 
спины, направленные противоположно. В обоих случаях получена зависимость того же 
самого типа: ὦ (k)—k?, что согласуется с макроскопическим поведением обоих родов MarHe- 
тиков. Формулу o (k)~k для антиферромагнетиков нельзя получить из выведенной за- 
висимости ὦ (К) для антиферримагнетиков посредством обыкновенного перехода к пределу 
5,-5;, т.е. рассматривая антиферромагнетизм, как особый случей антиферримагнетизма. 
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The separation of angular variables in tbe relativistic one-time equation for the 
two-fermion problem in quantum electrodynamics is carried out for the case of one- 
-quantum interaction. The momentum representation is used. А system of integral 
equations in one variable is obtained. This system reduces to four equations in the case 
of j = 0, and to eight equations for j Z 0. 


1. Introduction 


In a previous paper (Królikowski and Rzewuski 1955, referred to hereafter as I) 
the relativistic many-time equation of the two-fermion problem in quantum field theory 
was transformed into the one-time equation possessing the form of a Schrödinger 
equation. The operator playing the role of the potential in this equation is, in general, 

‚а complex quantity. The imaginary part of this "potential" is due to the instability 
of the two-fermion states and determines the life-time of these states. Using in the 
kernels of the original many-time equations the functonis 5" (as in Günther's paper 
(1952, 1954)), we obtain the single-particle theory of fermions; using instead the func- 


tions = p (аз in the Bethe-Salpeter paper (Salpeter and Bethe 1951, cf. also Schwinger 


1951)), we get the symmetrical hole theory. | 

The one-time equations of the many-body problem in quantum field theory were 
further discussed (Królikowski and Rzewuski 1956) from a more general point of view. 
The resulting equations describe bare particles; if, however, the renormalization 
procedure may be carried out (Matthews and Salam 1954, Johnson 1956), they become 
the equations of physical particles. In these considerations (Królikowski and Rzewuski 


(321) 
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1956) we have used the initial condition corresponding to the single-particle version 
of the theory of fermions. Using the initial condition of Feynman’s type, one gets the 
symmetrical hole version. 

The next step should be the solution of the equations and comparison with expe- 
riment. In the present paper we solve the angular part of the two-fermion problem and 
use the solution for the separation of angular variables. The resulting set of integral 
equations in one variable is then investigated further. 

In virtue of I (cf., in particular, Section 4. of paper I), we have the following 
equation for the two-fermion problem in quantum electrodynamics if we consider the 
one-quantum interaction only: 


[W — H, )- ihd,) — H, (— ТЭФ (zi, x9) 
= | dx ds x3 V (а, 77 (1,1) 
where 


Н, (p) = y? (icy;p + т; 2) = ca; p + В; m, c, (y) = ( iba, B), (1.3) 


V (z وند‎ χι, X2) 
ای‎ 1 1 = > ; c SE 2 2 > — ہے‎ => 
= 1 hc [yP AP (αι, xa; м1, zl + γ А (αι, αὐ; Xis 22)] Уз” Уф (1,3) 
and = E — ir. 
Using the kernels of the many-time equations in Günther's form, we get 


> 


V (ху, хь; αγ, αρ) 
--- e i dx, уф Sy" (αι — x4) в УР, (x, — ας) S, (x1 — 21) وگ‎ (xa — x) + 
+ f d νῷ SP ба — x) AP VP Dis — Ὁ S, (а — x) S (а; — 23 ] y у®, 
(1,4) 


where the upper sign corresponds to charges with opposite signs and the lower, to 
equal signs. Taking these kernels in the Bethe-Salpeter form, we obtain formula (1,4) 
in which S;* are replaced by iS,,. The functions D, = (1/2) D, and S}; = (1/2) Sp, 

۹ 1 1 


are the Feynman solutions -of the equations kia 


] D, (x) SCH δία) and (8-х). Sale = — i (x), 


and 


РОЗ, where 00 κα 


is the retarded solution of the equation . 


(2+ x) Sf" (a) δω. | 
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Making use of the formula 


2 z Hi(—ihô) t 


5, (α) = e i ô (x) y (1,5) 


we may rewrite (1.4) in the form 


= 2" = E 
V (%1, χο; Xy, Хо) 
0 


i - i E 
z H,(—ihd,)t = =) — Hi(—ifio)t 
Ee i е ГО, (αι — х, ct) е a eo 


—oo 


i ea i 2 
P Ν᾿ : ГР, (x = <a, ct) e 2 E |, (x1 — x4) Ó (Xa — x3) 


V (αι, хә) д (αι SE 24) б (xg = κο) )1,6(‏ کے 
where‏ 
D= γῷ γῷ γῷ γῷ = 1 — ааз. (1,7)‏ 


Taking the kernels in Bethe-Salpeter's form and using the formula 


hà : IUD 

δι) = A, (— #9) e (t) — S (x) where e (= |_1 iu 
A کال‎ #9) | + H, (— #9) {+ اپ‎ 
` на ι΄ 20 


we get instead of (1,6): 


V (ху, αρ; 34,509) = 
0 0 7 
1 => 1 > 
= z Н, (193) t - +H, (—ihd,)t 
= F іе?с ως ο OTE a 2 = 


—oo 


D, (αι — х, ct) e 


3 5 e = 

- 40 (—.ihd,) dt EI (iid) 2 гр, (б E ЕТ E HB )- t 
0 m | | 

+ A® (— ж) J dt “ΠΕ Гр, @ = ἄν ole Е На (- 01 
aie | | : 

E (—ihð,) t =н, T 1 : 


— AO (— ihà,) J dte ГР, (x, = προ) ο 
0 


x äis - ху) Ò ) - м) = VG ود‎ Ô (1 — 1) δα (16) 


конт: ems 
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2. Transition to the momentum. representation 


First, we perform the calculation for Günther's case. 
Equation (1,1), in the momentum representation, takes the form 


[7 — Н, (p) — Hy (РӘ) 6 (p.p) = | Фр: do pa V (Pis Pi pis ۶ CL 
where 
Ф (Py pj) = | dy αι ὧν κα < ٦ 
V (Py» Pas = | ἄνα den < Pi Pa | а 2 V Gs) 
< («ж [Рур >; (2,2) 


here we have 


St 


E 
< Z) ль | Di Da > = Dahl u, (pj uz (pa) е" 


(Н (P) w (D) = Е (р) u, (p), E;(p) = £ eV p? (то), uj (pu لچ‎ 


1 
— Di Xa 
ef ۱03 


Thus, because of (1,6) and (2,3), and using the formulae T 


178] : 
Du = 5 Gs [ dk ó, us ей, 


8.0) =< J dt oF im — 8 0 ш 


ο. | \ | z ex x ΣΣ 
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where 


i ὑπ 22 ох x 2 d = 
.V (p, p') = + = uj (р) u5(— p Tu (p) w (— р) x 
1 
"e )2h |p 4 Ç : 
л — р’ an, 1 , 
Р-Р (6 В) η (B = в) 


; (2,6) 


here we have used the notation E; = E, (p) andern, E; (p). 
We note that equation (2,5) in the position representation take the form 


[E — H, (— ihd)— H, (ih8)] Ф (= f & V G, z) @ G), (2,7) 
where 
Φῶ- J hp > ١م‎ < PY 
وھ‎ ۹9 9۷/٦٣٥٣۲ (2,8) 
here we have 
i++ 
> > ET => d 737027 
< х|р > = (zh) * u (p) us (— p) е (2,9) 
In the static approximation the kernel (2,8) has the wellknown form (in (2,6) we 
put E; — E; — 0): 
Kare! 
4л |£ | 


Equation (2.7) is then the Breit equation with instantaneous electromagnetic inter- 


Noe eene RE -γῷδῦῶ ο ΄ (2,10) 


action. 

The wave function ® (p) in equation (2,5) is labelled by the unit spinors un and 
ш: Ф (р) = (Φ, иг 
(Ф. ( р)) labelled by the spinor indices a, and ag, similar to Ф (x) = (Фл, (x)). 
и: the transformation 


Ф„„(р) = У ut) (p) u® (—p) Φ, (р) 


Uili: 


( p). We now go over to the new wave function, Ф( ϱ) = 


and using the formula 


Y: uf? (B uf? = م۸2‎ Gere (р) = ATE), رن‎ 


= 2 | E; (pH 
we obtain from (2,5), an equation of the same form as (2.5), where now 
ο e 1 1 1 (H,-FE)T Ge + Ei) 
V (p, p') = + iuc Pl 4 > : ; 5 К ar 
sm | E Ei (|p — p'| + — را‎ 1) 
"n (H, — E») ا‎ TE) - 
| (2,12) 


5 E (Ip - pee MEL в) 
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here we use the notation 


H,= ο. НС, НН) HH) 


and 


у, f Gs E) = 7(Е,, E) مر‎ (LE. E) - f (E, E) Ef (- Εν — E). 


Introducing the notation 


Te 
sehe se Ele л e H, (p), 


τ. ο 1.15 اج‎ ee He ee 
Е PEI 


we may rewrite the two-fermion equation with one-quantum interaction in the form 


> d. p? 
(e — h — h) 9 (p) =F, | = 
1 2 4 (2л)? |р — p’| 
sign | ἒι 81 (Ip = p| Sp £i — εἴ) 82 ez (Ip 0] کے و‎ 82) ds í 


where a ss ei /4лйс = „1/137“. 


In the Bethe-Salpeter case we may use a calculation similar to Günther's case. Then, 
instead of (2,4) and (2,6), we get, respectively 


RU TP DA n je? ar EIA [24 T T Se - T, Ké 
V (Py, P2; P1, P2) = + απ “1 (Ра) u (py) Г ui (pi) uż (p2) ó (py + py — pi — pa) X 
+оо 


1 ره‎ (р) و‎ № 
x — | do d 6 | B ape TOL eda 
m lo, (ρι)| + [0 pi c? PNE 


w, (ру) e э о» (pj) > Sp A? 
P (GH (w (ру) — @ (p1)) — ο] + ΡΈΕΙ Tm δ᾽. E η) = = ον] ὃν X 
is σης, e (o оз (ра) -- w2 (p3) -e| (2,4) 
and 


S EA, e? Km * = DR 4 ed 
V (p. p) = + az “1 (p) и (р) T ui (p^) uz (—p') х 


E 
1 1 22 
EE 1 
Оа) рр | . 1E S 
کھت‎ SC 1 میرف‎ 
|р р | T ο LE, (Ei Ei) 
Ey 
tmc ome Gei 
|р — p'| + — cz (E, — EN 
с |El 
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Formula (2,6") takes, in the static approximation in the position representation, the 
well-known form 
7 ο). - е? (1) “2H, (2) pope (1) дү AQ کر بھی ور‎ 
(x x ) جح‎ 4 [4% (—iA8) AY (19) — AZ’ (—ih8) Л (ih8)] e Ó (x — x’) 
х 


= V (x) ó (x — x). (2,107) 
Instead of (2,13), we now have the following two-fermion equation with one-quantum 
interaction: 
te E SE ىف‎ 9 ٤ 
ол) | р-р’ 
x δὴ (hy + ει) Г (4 + εἴ) 5 
sign o 2 - 2 à! 2 | 
5 Ene Ре k ер ا‎ 
&l 
٤ hs'-- 6ο) F (he -ἷ- ες - 
SL i ( 2 2) ( 2 SX 2) Ф (р) (2,13) 


; > pr, 0€; ^ 
£g €2 Г ашы щш» = 3 
2 


Both equations (2,194) and (2,13) may be written in the common form 


(e — hy — №) $ (p) 


a Al 4 ү ; ў 
πιο edes 8 y [(h, + ει) Г (hi + 81) M (e1, εἰ) + 
sign 
+ (ha + εν) D (k; + єў) M (e, e2] Ф (р), (2,14) 
where 
(1p — pi ЕЕ’ [lp es + e — 2] in Günther’s case 
M (e, =) = \ in the Bethe-Salpeter 


€ Ce =? , E E € τ 
2 |р —p е Lus. uin] dud (2,15) 


3. Solution of the angular problem, 


As a first step towards the solution of equation (2.14), we shall try to separate the 
angular variables. This is certainly possible due to the invariance of our dynamical 
system with respect to rotations in three-dimensional space. This invariance secures 
the conservation of the total angular momentum of the system, 

J; = = ih M, M, = m, + $, 
Го 9 و19‎ 2, 
кыт فور موا یسا تر ےھ ور دی‎ 
ШГ Ор سور تن‎ wk У) | 
@ = 1,2,3; i,j,k = cyclic permutation of 1, 2, 3) (3,1) 


&- 
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and, consequently, its commutability with the integral operator of equation (2,14). 
We may therefore consider a system of three commuting observables (there are three 
degrees of freedom left after separation of the centre of mass variables) consisting 
of the integral operator of equation (2,14) and of the differential operators 


Љ and Л} = Ji + Ji + /$. (3,2) 
The full system of equations now consists of (2,14) and of the two equations 

Ja ® (p) = hm % (p), (3,3) 
J? (p) = hj (j + 1) 9 (m. (3,4) 

We easily verify the commutability of the ιο (9.2) 
Lie 71 ο (3,5) 

which is a simple consequence of the well-known relations 
[M;, Mj] = — M, (3,6) 


It may be interesting to verify also explicitly the commutability of J, (and there- 
fore, also of J?) with the integral operator of equation (2,14). For this purpose, we 
first calculate the commutators of M; with h, and hg. From (1,2) and (3,1) we easily get 


[M;, hy] == [М,, Πιο] =0 (i a I, 2, 3) (3,7) 


Commutability of M; with an integral operator of the form f d, pr F( р, p) means that 

M, f dsp’ F(p, p) 6 (p) = f ар F (p.p) M; 0 (p), (3,8) 
where M; may be obtained from M; by replacing p; by p; and E: by = x Ta 
show that relation (3.8) is, in fact, true for the integral operator of equation (2.14), 
we first note that 


[M, Г] = 0 (i=1,2,3) (3,9) 
and, due to (2.15), 
.9M (s, =) ЭМ (e, ei , ЭМ (e, ғ) ‚ЭМ (e, ei 
Р Pn o pug اتیک‎ ng М δρ, کٹا‎ 


Because of (3,1), (3,7), (3,9) and (3,10), we may now write for the first part of the 
integral operator occuring in (2.14) 


M, f d p' (hy + ει) T (hy + ὁ) M (ει, εἰ) 
= Í da p' (h + εἰ) T (m; + 5) (hi, + εἰ) M (ει, εἰ) 
e if ds р’ (hy + e) Г E (hi + ει) M (в, εἰ) — 


ры (2M (S εἲ) 20-0 x 
- @ + & (uM ο 121 | 
ο a or (3,11) 
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The vector space under consideration consists of functions @ (p), which are 


quadratically integrable and therefore vanish for | p | > >0 (the integration domain 
is the entire three-dimensional space of the momentum variables). Thus we may 


carry out in (3,11) an integration by parts, the contributions from the infinite boundary 
being dropped. 


M; | ds p' (h + ει) T (hi + e) M (ει, εἰ) D (p) 
= [dp (h + ει) TM (ey, e) M; (i + 6) $ (p) 
= [dp (h + e) T (h, + εἰ) M (εν, εἰ) M; % (p). (8,12) 


An analogous relation holds for the second part of the integral operator in (2.14); 
this completes the proof of the commutability (cf. (3,8)) of this operator with M, 
and, therefore, with any function of M; as well, e.g with M?. 

Equation (2,14) is inconvenient for direct treatment insofar as it contains angular 
variables both under the integral sign and outside the integral by the intermediary 
of the operators A, and hy. We shall remove this dependence on A, and وم‎ by means 
of a unitary transformation 


@ = UO, F-UFu, (3,13) 
where F is an arbitrary operator and U is the unitary operator 
= ett? езе, $) (3,14) 
9 and ф being the angular coordinates of the vector p in an arbitrary polar coordinate 
system 
ру = p sin cos g, 
рә = p sin д sin g, (3,15) 
рз = p cos д. 


The set of equations (2,14) and (3,3—4) is transformed by (3,13— 14) into 
(e — 7 = hs) Ф (р) 


| 1 = NES ; 
= ш | Убе) uro S res) + 


+ (ig + ود‎ UT'U-v (hi, + ej) M (a, ell O (PI (3,16) 
J, 9 (p) = hm (p), (3,17) 
J? O (p = hj (j + Ὁ) 9 (p). gen 


*) Transformations of this type were used e.g. by Sommerfeld in the case of the hydrogen atom 
(A. Sommerfeld, Atombau und Spektrallinien, II. Band, 4. Kapitel and Zusatz 15) and by M. Günther 
(unpublished) in his treatment of the relativistic two-body problem. We are very much indebted to Dr. 
Günther for showing us his manuscript prior to publication and for valuable discussions. 


k. και, 
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In (3,16) D and hs as well as h, and kh, do not depend any longer on the variables 
9 and g. Let us verify this in the case of h,. The dependence on g and occurs here 
by the intermediary of the factor 


у p; = p [ sin ὃ (y? cos p + y? sin φ) + y® cos 9] 


1 (D ay 1 (1) 


-- yi $9 - 3D Pv 
κ ο ας р (y? sin ὃ + уу) cos 0) e 
) 1 TED 1 050 
RAM MUS رج‎ MUERE DU IU RO k 
2 ` e 3,19 
=e 6 py? е (3,19) 


The spinors y occuring in (3,14) commute with y®. The products aD) yf? as well 
as y yf 2) (i, j= 1,.2, 3) commute with. yf? and yQ. We have ШОС 


hı = U γῷ (ір, γω + т, ὃ Ότι 


2 © y 9 a γῷ fe το = نار‎ D y e zip γῷ) 9 — 
(3,20) 
and finally 
hy = y? (ру + то). (3,21) 
Similary 
ha = y (— ipp? + my). | (3,22) 


The commutation relations (3,7), (3,8) and (3,12) are, of course, not affected by 
the transformation (3,13). 


To find the dependence of the eigensolutions of (3,16) on the angular variables 
9 and g, we first have to solve equations (3,17 — 18). For this purpose we have to 


find the transformed operators CS and CH This is very easy in the case of ike Indeed, 
we have in polar variables 


M, = — + 53 (3,23) 


and 


(3,24) 


que 


In the case of J, the calculation is more involved. We have 


| д . д 
М = ο ιτ می‎ 


туз д | 
М» = — sing ctg 9 bj + cos p >» + $2. (3,25) 
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With help of the relation 
evi Yk* — cos y + y, y, sin у 


we find first 


д д 
Zem ۸ھ‎ 
д д 
-1 
TT a 


U sı, U-I = s, cos д cos p — s, sin y+ ss sin 9 cos 0, 
Us, U-I = s, cos à sin 9 + s, cos p + sz sin ὃ sin g. (3:21) 


Transforming now the formulae (3,25) we get, on account of (3,27), 


Y д д cos 

Μι = = N, — Sing a + E 58, 

— к 2 д зїп ф 

M, = — sin g ctg д 25 + cos ф БУЛ + TW (3,28) 


and, with help of (3,24), 
Mola. 2 1 22 д 5 
M; = And 28 cr a L dn? = — 2 sg cos РР + 53]. (3,29) 
Having thus obtained the transformed operators M; and М? ‚ we may proceed 
to the solution of equations (3,17 — 18). The first of these equations (3,17) takes, 
on account of (3,24) and (3,1), the form 


90 . - 
E" = im Ф (3,30) 
and possesses solutions of the form 
Ф (p, 9, p) = gu z (2, p). ; | (3,31) 
The boundary conditions Ф (quu py Ф (p + 2л, ὃ, р) are satisfied if m is an 
integer | 
nri, Е2.... (8,32) 


Introducing (3,31) into (3,18) we get, on account of (3,29), 


o д 1 : 
ET DEE Ee m? — 2 im э св д + [26,5 


Let us first solve equation (3,33) for m — 0. In this case if we write out s2 explicitly, 
it takes the form 


1 2 

quer 12 Ë Пр 7 
sin 9 20 (+ i 5] Male IEEE ای‎ (9, p) 
= -jG- 0x (8, р). (3,34) 


a. id, 
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Here use has been made of the fact that 


1 
vec ες vP) (3,35) 


is an idempotent operator 
1 
Ë (руи? И = σα- ym 2 ҮР уф). (3,36) 


Equation (3,34) possesses the form of the differential equation of spherical harmo- 


1 
nics with the upper index — s2 — 3 (1 — yP γῷ VO у!) and with the lower index j 


2 
Z? (cos 9) = P; 3 (cos 9). (3,37) 
The meaning of (3.37) becomes immediately clear is one considers that the two mutually 
1 " 
orthogonal operators 5 )1 - γῷ A γῷ ») апа > a + yQ у γῷ) у) are idem- 


potent and possess, therefore, the eigenvalues 0 ce 1. We may thus write instead of 
(3,37) 


РД (ооб) 2 (1 + У ٹر شر‎ у) PR (cos) + — 


E NT — γῷ yf P y y) Р} (cos 9) (3,38) | 
Due to the relation | 4 | 
c AUI ыш — 

e (EAR $92 =: cage SS: 

we may write (3.38) also in the form | SEE ge py es co m 


E It^ Eer z pne а 
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Due to the rotational invariance of our equations in three-dimensional momen- 
tum space, the equations for the radial component cannot depend on the value of the 
quantum number m and therefore we may already use the eigenfunctions Z of 
equation (3,34) for the separation of the variable 9. Nevertheless, it may be of some inte- 
rest to know also the eigensolutions of the general equation (3,33) with m an arbitrary 
integer. These eigensolutions were derived by M. Günther (cf. footnote on page 329) 
from Z? as a generating function by changing the pole of the polar coordinate system. 
The result is 


m 1 
Z; (cos 9) = Б (La: yo y у γί) zë 


1 1 , 
SETTE )1 — y? уз у уш) x 
и А 
x [55 سے‎ Hyd 25] р" (cos 9) (3,43) 


Thus for т = 0, we obtain the functions (3,40). 

It may be verified directly that the 2” are eigensolutions of equation (3,33). 
For this purpose it is convenient to write the second order differential operator occur- 
ring in (3,33) as the square of a first order differential operator 


E 9 1 A | 
ET т - (3,4 
sin ð 8 Є d 6 25] sin? д S SR и (3,44) 


where 
= Dy od 3,45 
ou (5 7 gi due SE |. | (49 
It is now easy for us to verify that اھ‎ is an eigenfunction of K. Due to the orthogonality 


1 
of a +P YP yP УФ) and — ss= (1 - y Y? YO AEN, we get 


т 
КА" = TIE )1 + Y? y? vi ^X (5 o EM » 90 | 


Uer 1 ) m? ke 
+ ба جا‎ yU у y? y$) Е Rt sare) | 7 
=үў+1) Fry? γῦ y? (3,46) 
and therefore hy Pea 
E he у (3,47) 


4, Separation of the angular variables 


Having thus obtained the eigensolutions of the angular problem 6 17 — 18) 
we may now proceed to the separation of the angular variables from the integral ` 
equation (3,16). As already stated above the result of the separation cannot 


depend on m, and therefore we may use for our purposes the functions 2} (соз 9). 


554 W. Krolikowski J. Rzewuski 


We write the solution in the form 


Ф (p) = Z (cos 9) f (p); (4,1) 
we introduce this into the integral equation (3.16), multiply on the left with Z? (cos 9) 
and integrate over the solid angle (dw = sin 9 d dg). Due to the commutability of Z? 
with 7 and hos we may write the result in the form 


1 πο mn πω, 


со 


== = fær = 201 رآ‎ E ED PK; ER 
0- sign 
25 (h, + εν) ГК; (es, εὖ) (h2 + =2)] f (p?) (4,2) 
with 
K, (e, ε) = i! do dear Z? UM (e, ei U'-1 ZW. (4,3) 


The functions Z? form, according to their definition (3,38) and according to 
(3,41), an orthogonal set 


J 49 зп ὃ Z2 ZR =0 forj عو‎ k (4,4) 
0 
For j > l they are normalized to unity 
f 4950 7972 =1 (4,5) 
0 
and for j — 0 
А 1 
I sino 20 Д0 = 7 (1 + γῷ γῷ ye у). (4,6) 


0 
We may therefore replace, for j کو‎ 0, the integral on the lefthand side of equation 
(4,2) by 2л. For j = 0 we get 2л times the factor (4,6). This factor, however, occurs 
on both sides of the equation in connexion with the appearance of Zo 0 = E (ee: 


+ yD yP 40 y) Ро Tt commutes with h, and SS and is idempotent. It may Im 


be joined with f (р). Thus, for j = 0 as well as for j Z 0, we may write equation (4,2) 
in the final form 


(e — hy — وط‎ f (p) 


со 


эе ls NC + а) ГК, (ει, εἰ) (hy + εἰ) + 


sign 


Be (ha + ë) T' K; (ез, εν) (h, + e)l f (p). (4,7) 
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The angular variables no longer occur in this set of equations. It determines the 
energy eigenvalues of the relativistic two-body problem in the desired approximation 
(one-quantum interaction). 


| 1 
The fact that for j = 0 the projection operator > (1 + P γῷ γῷ) у?) occurs 
аз a factor in the solution f (p) shows that the lowest state of the total angular momentum 


may be realized only in such a way that the orbital momentum as well as the total spin 


(1) 0) „,(2) 


13 zero ib (yi yi Yi 


у) f (p) = d ‚ he. the particles possess opposite spins. 


For further treatment of equation (4.7), we-need, first of all, the explicit form of 
К, (е, el, Instead of К, (е, ε) we shall calculate the expressions 


Ky (e, ε) = | dodo! Z? UM (e, ғ) U' ZY (4,8) 
and we shall show that 
Ky = رآ وڈ‎ (49) 
This result will at once confirm our statement that the functions Z? bring equation 
(3,16) into the form (4,7). In fact if we introduce (4,1) into (3,16), multiply on the left 
with Z? (L Z j) instead of 7? and integrate over the solid angle, we get zero on the left 
due to (4,4), and zero on thé right due to (4,9). 


To calculate (4.8) we first note that М (e, =’) depends implicitly on the angular 
variables through the intermediary of 


|p — р = V p? + р? — 2pp'cos 0 (4,10) 


(cf. (2.15)), where cos 0 = cos д cos δ΄ + sin 9 sin 0" cos (p — g’) is the angle between 


p and р’. We may therefore expand М (e, =’) in a series of Legendre's polynomials 


M (е, г’) = 25, М, (е, =’) P, (cos 0), 
k—0 


M, (e, &) = [ M (е, =) P, (cos 0) sin 0 d 0. (4,11) 
0 
Inserting (4.11) into (48) we obtain 
Kj УМ) 7ھ‎ 00 д 6 
Š Bere ¢ 4H | 8 
A= f do do! Z° ОР, (cos 0) U'1 Z”. | (4,13) 


Using finally the summation theorem of spherical harmonics 


` 1 m= +k 
P, (cos 0) =- 


жи; 
Е, 


РЇ (cos 8) PI (cos ϑ') e" 0-9 (4,14) 


bf, A ech 
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(Р;”" = (— 1)” Pf), we may write (4,13) in the form 


iy 


A 
vm / | СЯ κ Lee; 


We first carry out the integration over the angles ф and φ’. On account of (3,14) we 
have to evaluate the following integral 


0 m ei" (9 — 9) рт Tr! —1 
> ees UP? PRU ZI. 419) 


л ہو‎ 2л 
fa J do’ eim 9-9) es 0-0) = zu dg e"? e»? (4,16) 


With help of formula (3,26), or expanding the exponential and using the fact that 
— 52 is an idempotent operator b (3,35 — 36)), one easily verifies the relation 


2 > DÉI ار‎ D 4 42 42) φ 
e = 
Je E 1 2 γί (1) γω уб 2). y2) wen =q UM yb y γί (2) y) cos و‎ + 
E 4 
Ed 0 πο. yt y) sin φ. | г) 


BE (4,17) into (4,16) and carrying out the жыш one zeit 
j dp Г dq’ gne ος ЕСЕ AE y, y y y) Ons 


1 
4 2 (1— у E у? E 5 On Z^ x к kt Ab ۱ 


Lui Y + y y) — (ó, 
Së Kate, "x a 
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We are thus left with integrals of the type 
I, = | do sin 9-20 et” 20. (4,20) 
0 
Using (8,1), (3,35 — 36) and the mutual orthogonality of the operators 


(1) (2) „,(2) 


ο Пу YP y? v6 


), we obtain 
I 2 dis | 5 I+ yp yp تار‎ y?) P? P 4- σα — γᾶ) y γῷ) Pj "| x 
0 
ДЕ (939 9 30) у ü — yD P γ᾽ УР) οὐ JE 


` 1 
«f dh οἱ in; (yD γῷ 9 γῷ) PPL 0 — 777 92) oi x 


1 ۱ 
x FOP AP + y yË) sind. (4.21) 
The integrals may be evaluated with help of the well-known recurrence relations for 
spherical harmonics. In our normalization these relations take the form 
LEL Аа СС ИРА 
V Gj +1) Gj + 3) a = Т (27 + DE 


cos ὃ P; == 090007 + КОЕ вр: Рем. 
(2j + 1) (3j + 3) (2j — 1) Qj + 1) 


sin OP? — — G+)U+2) p Pla Vai GS DZ „pl, (492) 


cos 8 P$ = Pia ct 


(27 + 1) (2j + 3) (2j — 1) @j + 1) 
Putting (4,22) into (4,21) and employing the orthonormality relations (3,41), one gets 
after some calculations 


i Ki ا‎ == y® yí (1) ye y) д 4 


a σοι (НУР 9 yy) [урт у تاپ‎ VG 0 9 2] 
l0. 499 0 2) [VG 12) + у 3? VjG E 3] | ME + 
n E (1—y γῷ уб تا‎ (1 + γῷ ناپ‎ γῷ у2) [j ART) γῷ Wr - SE + 


+P PNG DFF 


1,k 
rp VOF ЕЕ pP 
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From the above one obtains after a rather lengthy, but straightforward calculation 


1 Е 
I ly = 5 (1 + ak? ah? ys 7ο) да ды + 
1 ET (2) 702) TE. 
+ а PAP a + BP ۷٣٣ 
(2) 


1 d ail yP у? le (1). (0) 2 E ó S: 4.94 

БЕ 5 Hoc EDUC pem ИЗ, Vid Е 1) k,1—1 ἮΙ + ( , ) 

With help of (4,24) we may easily calculate A, „ (cf. (4,19)) and introduce the result 

into (4,12). Due to the occurrence of the Kronecker symbols, we obtain, finally, 
relation (4,9) with 


K. 1 М, (е, е) 
air. 0r کر‎ ив) mE, 
)27( 2 us 
2 
1 1 D νῷ) 4,02 4,(2) ES М; + 1(в,=’) 
FSD) [1-285772 85 E von) 0 
2 2 9 -E d rod 
2 
1 у ys ym nu (1) (1) — Mie €) 
= AT LE Xd esi μμ A Mall as y 2 کک ن‎ 
SE E ا رت‎ 77٣ . 1). (4,25) 
oF 


The functions M, (e, г’) occurring in (4,25) are given by the formulae (4,11) and 
(2,15). The integrals may, in principle, be evaluated by elementary integrations, but 
we did not succeed in finding the general formula for arbitrary j. We shall write out 
these integrals explicitly using Günther's form of interaction (cf. (2,15)) and evaluate 
them for the two simplest cases: ; = 0 and VENE | 


M. : 1 P, (cos 0) sin 0 š 
k (68) = — MA UAR oO N (0,06) 
ё J pš р’ — 2pp' cos 0 + (e — гүр? + р” — 2pp' cos 0 
0 ‘ 
/l1 1 pt p+ ee 
MaE e) у с in | 
o (= е) 2 εε pp’ In ск 


; Я & — ε’ 
men- i-i 5 is 


Ера (ee) Ae pes wer ] 


; - ; 4,27 
2 pp р-р Γε-ε (120 


9. Reduction of the system 


The integral equation (4,7) represents a system of sixteen equations for the sixteen 


components of the spinor f (p) (cf. the definition of o ut p) on page 325). This system 
may be reduced, for j > 1, to eight equations for eight unknown functions, and for 


= ہم 
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J = 0, to four equations with four unknown functions. This is due 10 the fact that 
for j > 1 the spinors y? and y® occur only through the intermediary of their product 
yP γῷ) (cf. (4,7) and (4,25)). For j = 0 the terms in (4,25) containing у y{ vanish 
and the resulting equations contain another pair of spinors occurring only in the 
form of their product, namely, γί) and у. 

Thus, in the general case j > 1, there occur in our equations (4.7) the spinors 


Ио, уф, yP, P, Pp, бл) 
satisfying the following commutation relations. 
bé yo], Se? 2δ,, 
[o У = 2Ó uy (for и, dee 1, 3, 4) (5,2) 
[γ9) yP, XP], = 0 
[ν9» УР, УУЛ. = 0 (for » = 1,8, 4) 
(o, b] + = ab + ba). (5,3) 
These relations are satisfied if we represent the spinors (5,1) as products of the second 
rank Pauli matrices ; 
1, of, 05, σᾶ; [o?, of], = 205.0? оў = io} ; (5.4) 
corresponding to three different two-dimensional spinor spaces. To distinguish be- 
tween these spaces we denote the corresponding Pauli matrices (5,4) by σι, о; and 2, 


respectively. 
We now put 


y = gs 0, = o$ x Of 1 = (o8, 01,5 д). 

y = o, 0, = 05 X of X 1 = (of, و‎ δη), 

yP = 2,0 = of X of X 1 = (07, 01,4 д.) 

yP = 0, djs = 1 X Of X 08 = (б 01, 93.) 

y = 01 У = 1 Χ of X оф = (0,5 07, 05,,). 

ak) дү 225 = 1 x of X 0 SZ (0,5 013 с 

γῷ у? = 05 = 1 x of X 1 55 (0,5 о, ga, 2 | (5,5) 
One easily verifies, with help of (5.4), that the spinors defined by (5,5) satisfy relations 
(5,3). ا‎ 

We may now express all spinors occurring in our equations by means of the products 

of the second rank matrices oy, 01 and .رھ‎ The function f (p) appears now as а spinor 
possessing three indices В, ὃ, η corresponding to the three two-dimensional spinor 
spaces. Each of the indices takes on the values 1 and 2. Equations (4,7) represent, 
therefore, a set of eight equations for the eight components fgg, (p) (В, 9, n = 1, 2) 


of the spinor f (p). 
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In the case j — 0 there occur only the spinors 
2): 503502. 4,0) 
Y, 9. 99. γῷ, AP YD. YE? YE (5,6) 
satisfying the commutation relations 


Ee ye by = 288 


bh m2, (μι بد ک‎ m 
bP PP Pl —0 (uv 34) б®) 
[γῳ γῷ, y], = 0 

ys @) اتا‎ =0 vue + 8680 (5.9) 


To reduce the system of equations we need now only two two- dimensional spinor 
spaces with the Pauli matrices o; and X; Putting 


ox = σῇ X 1 (of, д)‏ = لاپ 

оа = Lo ap da‏ نز 

у = 2 = 1 X оў = (δι 05,9)» 

rP = 22 = 1 X of = (дш 07,6) 

yy) = уру = оз 25 = o$ 03 = (08.8 08у) (5,10) 


we satisfy relations (5,7 — 9) and express at once all spinors in our equations in terms 
of the second rank matrices o; and Z;. We thus obtain a system of four equations for 
the four components fx, (р) (В, ó = 1,2) of the spinor f (р). 


6. Discussion 


We did not succeed, so far, in solving equation (4,7) exactly. This seems to be 
very difficult, if at all possible. However, equation (4,7) offers already certain advan- 
tages as compared with the original equation (1,1) or (2,14). Due to the fact that it is 
an integral equation with one independent variable only, one has better possibilities 
of applying variational numerical methods here than with the original equations in 
three independent variables. In general, the fact that one obtains equations with one 
independent variable seems to be the main advantage of the one-time formulation over 
the many-time formulation of Bethe-Salpeter and Günther, which even after separation 
of the angular variables (cf. Günther's paper, footnote, page 329) remain equations with 
two independent variables. 

Another advantage seems to be the possibility of applying stationary perturbation 
methods. E. g. in equation (4,7) we may use as the unperturbed equation the static 
approximation of (4,7) in Günther's case (cf. (2,10)). Putting in К; -۔ی‎ ё =0 wege 


со 


(eh Mm: و‎ i ἀρρ KE (p, p") др) 4&0 3 
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with ἘΣ (р, р’) given by the formula (4,25), where for М; one has to take expres- 
sion (4,26) multiplied by ee’ and with = — #' = 0 


со 


ae ^ Í Pr (cos 0) sin 0 d 0 
М, (p, p’) XX E J we = 2pp' cos 0 (6,2) 
0 


If it is possible to find exact eigensolutions and eigenvalues of (6.1), one may 
perform the conventional perturbation calculation to find the corrections to the eigen- 
values and eigenfunctions caused by retardation and, in this way, estimate the appli- 
cability of the static approximation (6.1). We hope to return to this question elsewhere. 
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This paper describes a two-sided spark counter for detecting thermal neutrons. 
The counter has cathode area of 7.5 x 11 cm? and is covered on both sides by a coating 
of boron. The total length of the 0.08 mm diam. nickeline wires which form the anode of 
the counter is about 3 mm. From the accompanying characteristics drawn for different 
concentrations of CO, in the air forming the atmosphere of the counter, it can be seen that 
as the CO, concentration increases, the counter efficiency increases, the plateau length 
decreases and the slope of the characteristic increases. The conditions necessary for the 
stable operation of the counter are given. 


1. Introduction. 


In recent years a new type of counter, the spark counter has been developed. 
This counter is particularly suitable for recording strongly ionizing particles. Since 
the first prototype of this counter was built by Chang апа Rosenblum (1945), some 
constructional changes have been introduced. From these changen: three basic types 
of spark counters have evolved. i 

1) A counter consisting of a cathode in the fofi of a plate and an anode in the 
form of one or more thin wires running parallel to the cathode at exactly the same 
distance from. it. 

2) A counter consisting of two еца plates, опе of which forms the anode and 
the other, the cathode. Е 

3) A counter consisting of a cylindrical cathode and a wire running along the ` 


` axis of the cylinder. This type of spark counter has a design similar to G—M counters, 


but works at a considerably higher voltage and is characterized by a different, type of 
discharge (Bella et al. 1948, 1949; Connor 1949, 1951, 1952; Eichholz 1952; Madansky 
et al. 1948, 1949; Müller at al 1953) ^ 

. Hereafter we will restrict ourselves to a description of the properties of the first- 
of the above. three types of spark. counters. Séi 


(343) 


344 J. A. Janik, A. Szkatula 


As already mentioned, this counter consists of a cathode in the form of a highly 
polished plate and an anode in the form of one or more wires of about 0.1 mm diameter 
stretched parallel to the cathode at a distance of 1 to 2 mm from it. Due to the applied 
high voltage (of the order of several kV) a corona discharge takes place. If an alpha 
particle passes through the area of this corona discharge between a given wire and 
the cathode (in a direction more or less perpendicular to the cathode), a spark occurs. 

The counter thus described is completely stable in dry air at atmospheric pressure. 
It was established (Payne 1949) that the addition of a small amount of carbon dioxide 
increases the efficiency of the alpha particle counting and even the ability to record 
protons and beta particles. 

This type of spark counter was employed by Savel (1952a, 1952b, 1954) to record 
slow neutrons. This application is based on the coating of the counter cathode with 
a layer of amorphous boron. Czownicka and Janik (1955) designed a counter similar 
to Savel's. All these counters had a relatively small cathode area (about 20 cm?) and 
had a number of wires, whose total length did not exceed 80 cm, stretched above 
one of the surfaces of the cathode and running parallel to it. 

In this paper, a two-sided spark counter for recording slow neutrons will be 
described. This counter is distinguished by its relatively large cathode surface (about 
100 cm?) and total length of wire of about 330 cm. 


2. Spark Counter Design 


The design of the spark counter for recording slow neutrons is shown in Figs. 
1 and 2. The basic features of its design are as follows: The cathode is a well polished 
` aluminium plate having an area of 7.5 x 11 cm? and thickness of 5 mm and coated 
on both sides by a layer of amorphous boron ὃ mg/cm? thick. For coating the cathode, 
a proper amount of boron was weighed out and a suspension of the boron in alcohol 
"was prepared. А small amount (1 to 2 drops) of Canadian balsam was added to this. 
Then the suspension was smeared over the entire surface of the cathode with the help 
of a brush, an effort being made to obtain as even a coating of boron as possible. 
After drying, a relatively strongly adhering boron layer resulted. However, because 
of the dielectric properties of Canadian balsam, the preparation of the layer in this 
manner made it difficult to obtain a corona which was evenly distributed along the 
wires, and therefore, the boron surface was then coated (by spraying). with a 0.1 u 
layer of metallic aluminium. 

The anode of the counter was constructed from resistance wire of 0.08 mm dia. 
wound between two brass posts fixed in such a manner that the cathode was located 
in the counter between the two rows of wires. Thus a set of 22 wires faced both sides 
of the cathode surface; a two-sided counter was therefore obtained. 

Plexiglass was used as the insulating material. 


The counter was placed in a metal box fitted with two glass plates whose surface 
was parallel to the surface of the counter cathode. The glass plates were set on rubber 


Fig: 2. View of the spark counter. 


i) screws for fastening the counter inside the gas-tight box described in the paper. 


the rigidity of the counters g) high-voltage input to anode A) screws for holding place the plexiglass slats 


e) plexiglass washers f ) plexiglass slats for increasing 


Fig. 1. Diagram of the spark counter: a) cathode, b) anode, c) brass frame on which the anode wires 


are wound, d) screws for holding the wires in place, 
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the counter was placed was supplied with two metal cocks 
with different gases. 


could be filled 


through which the box 


fastened to the high-voltage terminals by suitably chosen screws. The box in which 


seals and high-voltage terminals were mounted on one of them. The counter was 
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The following features of the counter should be noted, since they are decisive in 
determining the stability of operation: 

a) The wires must be closely maintained parallel to the cathode surface; the 
wire tension must be chosen so as to prevent deflection of the wire under the influence 
of the electric field. b) The counter must be kept dust-free, and the boron grains 


300 


200 


150 


100 


50 


Fig; 3. Characteristics of the spark counter for various concentrations of CO, in air: Curve I — 
air under a pressure of 1 atm. with 0% CO,; Curve II — air under а pressure of 1 atm. with 2% COs; 
Curve III — air under a pressure of 1 atm. with 3% СО»; Curve IV — air under a pressure of 1 atm 
SS S ' with 5% CO, | j | А 
cleaned before spraying, etc. c) It is recommended that the counter be filled with 
air (or air mixed with .CO,) from which the moisture has been carefully removed. 
Traces of water vapour can prevent the operation of the counter. A well-prepared 
counter has a background of about 1. discharge per 6 minutes. Any inhomogeneity 
of the corona glow (strong glows appearing at points along one of the wires, or different ` | 


> ۹ 
~~» 
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intensity of glow for individual wires) can lead to a considerable increase of background 
or even make the functioning of the counter impossible. If the above conditions are 
fulfilled, the corona will have a uniform intensity of glow on all wires. This appears 
to be an indispensable condition for the stable operation of the counter. Other con- 
structional features of the counter shown in Fig. 1, but not discussed above, are designed 
to prevent surface discharges, facilitate the winding of the wires, and so on. 

Fig. 3 shows the characteristics of the counter for different concentrations of CO, 
in the air filling the box of the counter. As can be seen, an increase in the concentration 
of CO, causes an increase in the efficiency of counting slow neutrons and, at the same 
time, shortens the plateau of the characteristic and increases its slope. 

The counter registers 0.025% of the neutrons falling on the surface of the cathode. 

This counter was used by the authors for determining the moisture of walls by 
using neutrons (Janik, Szkatula 1955). 


3. Electronic Equipment used with the Counter 


In Fig. 4 is shown the schematic diagram of the electronic equipment used with 
the counter. 


Fig. 5 is a photograph of the arrangement. 


ECL ft 6H6 А72! 
жеты 


E 

Fig. 4. Schematic diagram of electronic equipment: r, = 1M, το = rs = га = rs = 2.5М, r, = 10M» 

Fee 600 Q, Lá = 1M, سے وہ‎ 0.2 M, ho — 2M, Tu = 240K, г» = 500К, Cı = C, = G, = C, = 01 F, 
, С; = 100 z F, C, = С, = 16 pF, C, = 50000 pF, C, = 60 pF, Cu = luF. 
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Fig. 5. View of spark counter together with electronic equipment. 
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LETTERS ТО THE EDITOR 


SLOW NEUTRON DETECTION IN PRESENCE OF GAMMA RADIATION 


By Тлреџѕ2 0011۸07 
Institute of Nuclear Research of the Polish Academy of Sciences, Warsaw 
( Received June 5, 1956) 


I present а gamma-insensitive detector of slow neutrons consisting of: scintillator, 
931-۸ type photomultiplier tube, amplifier and discriminator. Although this type 
of photomultiplier is disadvantageous in scintillation counting it had to be used for 
reasons of availability. 

The scintillator consists of ZnS(Ag) used in form of a thin layer deposited on the 
glass bulb of the tube. The ZnS(Ag) layer is covered with a natural boron layer. The 
luminescence time of this luminophor amounts to a few microseconds. 

The scintillations engendered by alpha particles from B 10 produce large electrical 
pulses on the multiplier tube output. Unfortunately, pulses produced in the tube 
itself by a strong gamma radiation are of the same order of magnitude. Gamma rays 
do not produce appreciable scintillations in the ZnS(Ag) layer but they expell electrons 
from the photocathode and the dynodes. 

There is, however, a fundamental difference between both kinds of electrical 
pulses as the pulse duration is large for scintillation pulses and short for gamma 
pulses. In order to separate these two kinds of pulses a Miller integrator circuit 
followed by a simple amplitude discriminator was introduced in the discriminator 
stage. In this way gamma-ray background can be practically eliminated. The efficiency 
of this detector is close to 1%. 
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ROLE OF TRANSVERSE MAGNETIC FIELDS IN THE 
ELECTRODELESS GAS DISCHARGE AND IN THE JOSHI EFFECT 


Ву S. R. MoHANTY AND T. В. ΒΗΑΤ 
Physico-Chemical Laboratories Banaras Hindu University INDIA 


(Received June 12, 1956) 


The Joshi effect (Joshi 1943, 1945), viz., the instantaneous and reversible photo- 
variation (mainly diminution — Ai and, under definite experimental conditions, 
enhancement + Ai) of the current i of the gas discharge, is of fundamental importance 
for theories of the behaviour of electrically excited media and the influence thereon 
of irradiation. It is observed in numerous gases and vapours, elemental and compound, 
over a wide range of the operative conditions, e.g., the exciting potential V (and the 
attendent electric field E) (Mohanty, Jayaraman and Krishna Rao 1954), the gas 
pressure (Mohanty and Kamath 1948), etc., and is apparently a general property of the 


Fig. 1. Role of Transverse Magnetic Fields in the Electrodeless Gas Discharge 


discharge phenomenon (Mohanty 1955). A transverse magnetic field H modifies elec- 
tronic/ionic paths under E (Thomson and Thomson 1933, Townsend 1947). It was 
of interest therefore, to investigate in some detail the role of H in the electrodeless gas 
discharge and in the production of + Ai. 
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The electrical discharge was produced in a chemically cleaned soda glass tube 
(outer diam., 7.27 mm; inner diam., 5.57 mm) filled with pure oxygen at 26 mm Hg 
(26°С). Two external rings of untarnished copper foil (each 0.5 mm wide) separated 
through 2.8 cm constituted the electrodes. The gas was excited over 2—4 kV (rms) 
of 50 c/s frequency. The effect + A i (= ij, — ip, where i, and ij are the values of i 
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Fig. 2. Variation of the Discharge Current and the Joshi Effect with the Transverse Magnetic Field 
under light and in dark respectively) was observed on a crystal (1N34, Sylvania) — 
galvanometer (G) system (Fig. 1) with light (3700 — 7800A) from a 300 watt, 200 volt 
tungsten filament lamp placed at a distance of 26 cm from the discharge tube; heat 
radiations from the lamp were cut off with an appropriate filter. Magnetic fields H 
(0—750 gauss) of desired strengths were obtained by passing suitable currents through 
a solenoid consisting of 600 turns of enamelled copper wire (B. S. G. 24) per cm on 
a uniform glass tube (diam., 4 cm). The solenoid was branched into two equal parts, 
each on either side of the discharge tube, such that the produced H was at right angles 
to the normal electronic/ionic paths under E. : 

It is seen from Fig. 2 that i decreases progressively with H. Thus e.g., ipat 3.6 kV 
. decreased from 104 to 62 for an increase in H from zero to 750 gauss. The threshold 
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potential Г of the self-maintained discharge, characterized by a sudden increase 
in ¿ (Fig. 3), on the other hand, is increased linearly (inset, Fig. 3), from 2.89 kV 
at H — 0 to 3.13 kV at 750 gauss. 

The positive effect + Ai is small at low applied V and increases initially ceteris 
paribus with V (Fig. 4) to a maximum at V - Above یی‎ + 4 
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Fig. 3. The Do akdown Threshold Potential as а Function of the Transverse Magnetic Field 


diminishes rapidly and changes sign (to the negative effect — Ai) at an inversion 
potential V7. Beyond V/, — Ai increases with V in the range investigated. Whilst 
+ Ai at И < V4.4; max diminishes with H, that in the region V, A max — 2 
increases (Fig. 2). Furthermore, the maximum + Ai is smaller the higher the Н 
(Fig. 4). Thus, V, max and Г! values for the tube (Н = 0) were 2.75 and 3.11 
kV respectively. At 2.53 kV, + Ai was 13 for zero, and diminished to 4 for 750 
gauss. At 3.07 kV, on the other hand, + Ai values were.4, 6 and 10 for H = 0, 450 
and 750 gauss respectively. The negative effect — À i at V > VÎ diminishes with, 
and inverts at higher, H. At higher V, — 41 diminishes progressively, e.g., at 3.6 kV, 
from 13 to 5 over 0 — 750 gauss. The potentials V, A ma, and V7 are enhanced 
under H (Fig. 4), the former from 2.75 to 2.85 kV and the latter from 3.11 to 3.44 kV, 
in the range investigated. 

For occurrence of + Ai, Joshi (1946, 1947a, 1947b) postulates that (1) an ad- 
sorption-like electrode layer consisting of electrons, ions and excited particles is formed 
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under discharge; (2) light releases electrons from this layer; (3) пеш ion form- 
ation due to capture of the photoelectrons by excited neutral particles reduces $ to 
give — Ai, mainly as a space charge effect; and (4) uncaptured photoelectrons, and 
their secondaries produced under E if appropriate, cause + Ai. The effect of H is to 
produce a drift of electrons/ions at right angles to both H and Ё. Consequently, the 


20 
PO, : 26 mm (26°C) 
1N34 50^ 
15 
10 3 
1/2/3/4/25/6, N 
us SE 
η 0 
= Φὴ 
š 
10 τ 
° x 
| = 
О 
15 H (gauss) 
200 400 600 ` 
2 2l 3 352) 


Applied potential V (kV,rms ) 


Fig. 4. Influence of Transverse Magnetic Fields on the Potential Dependence and Inversion of the 
Joshi Effect 


removal of electrons by diffusion and wall recombination is accentuated (Loeb 1939). 


The observed diminution in i, as also in + Ai , with increase in H follows. In conse- 
quence of the electron/ion drift, an electric field Е” is set up opposite to the original E. 
The observed diminution of + Ai below V +4i шах › ЇЗ enhancement in the range 
У +4 max — Vj and the diminution of — Ai follow, as also the augmentation 
in V, and VI. | 


Grateful thanks of the authors are dus to Professor S.S. Joshi for his kind interest 
in the work. 
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ERRATA 


Acta Physica Polonica, 14, 135 (1954): 
Disintegration of the ®Be Nucleus in Coulomb Field by J. Sawicki 


Recently W. Czyz (Phys. Rev., 102, 1185 (1956)) has noticed that Guth and Mullin (Phys. Rev., 76, 
254 (1949)) erroneously omitted the interference terms coming from the different P — S and P — D 
transitions in the ?Be (y, n) angular distribution. Here the interference terms in the angular distribution 
were also omitted after Guth and Mullin. 
On including these terms the following corrections should be made: Eq. (13), p. 138, should read 
12 


| I, (P > D) |? fep (2,) = ов 1 4-3cos? y) | Аьр |? (13) 


redi 
Eq. (17), p. 140, should read 
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32 
Gh = 7/2 $20? 729 [ Rps J2 + | Rpp la + 2 Rps Rpp cos (9, —5 δρ) ln T' + 


1 1 
4+3 i Крр |? — Rps Rpp cos (0, m >) E (3 cos? 0,— 1) (1 — 72) + sin? бш r) | (17) 


Since the interference terms vanish after integration-over Ω,, the cross-section == do/dE, (Eq. (18)) 
remains unchanged. 


Acta Physica Polonica, 15, 49 (1956): 


Generalization of the Method of Supplementary Variables to Systems Composed of two Kinds of Particles 
by Z. Galasiewicz 


4ne? N, 4ле? N, 
On page 55 line 5, the expresion should read = 
47727262, 471726277, 
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